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Introduction 


This report summarizes research done under NASA Grant NAGW-2541 through September 
30, 1997. The research performed under this grant includes development and maintenance 
of scientific software for the GOME retrieval algorithms, consultation on operational software 
development for GOME, sensitivity and instrument studies to define GOME and SCIAMACHY 
instruments, consultation on optical and detector issues for both GOME and SCIAMACHY, 
consultation and development for SCIAMACHY near-real-time (NRT) and off-line (OL) data 
products, and development of infrared line-by-line atmospheric modeling and retrieval capability 
for SCIAMACHY. The European Space Agency selected the SAO to participate in GOME 
validation and science studies, part of the overall ERS AO. This provided access to all GOME 
data; The SAO activities that are carried out as a result of selection by ESA were funded by the 
present grant. The Global Ozone Monitoring Experiment was successfully launched on the ERS- 
2 satellite on April 20, 1995, and remains working in normal fashion. SCIAMACHY is currently 
scheduled for launch in early 2000. The first two European ozone monitoring instruments (OMI), 
to fly on the Metop series of operational meteorological satellites being planned by Eumetsat, 
have been selected to be GOME-type instruments (the first, in fact, will be the refurbished 
GOME flight spare). K. Chance is the U.S. member of the OMI Users Advisory Group. 



SAO performed the initial sensitivity studies for GOME and SCIAMACHY that determined 
their capability to measure gases including: 

• From GOME - O3 (including profiles and tropospheric ozone from the GOME nadir mea- 
surements), O2, O2-O2 (the collision complex), NO2, H2CO, H2O, SO2, NO3, CIO, OCIO, 
and BrO (global distributions). 

• From SCIAMACHY - the species measured by GOME plus, at least, CO, CO2, N2O, and 
CH4. Additional profiles will be available from limb and occultation viewing geometries. 


All of the measurements we proposed for GOME based on the SAO sensitivity studies have 
now been measured except for NO3 (due to the limited amount of eclipse mode data). The 
CIO measurements are currently at a rudimentary state, but will be improved substantially as 
radiative transfer studies improve (detectability of CIO under ozone hole conditions has been 
demonstrated and the instrument signal- to- noise is adequate for quite precise measurements). 
Thus we are quite confident that SCIAMACHY will provide the measurements as we have 
predicted, including the important tropospheric species CO, N2O, and CH4. 

The results of the various radiative transfer-related studies funded by this grant, coupled with 
development of advanced fitting procedures, currently permit SAO to make the most precise and 
accurate determinations yet of trace gas concentrations from GOME data. Figure 1 shows an 
example of a single BrO measurement from GOME, coincident with the SAO FIRS-2 balloon- 
based measurements (which include HBr and HOBr, as well as some 20 additional atmospheric 
gases). Such measurements are determinable from all GOME spectra, some 30,000 per day since 
mid -1995 through the present. Figure 2 shows an example of BrO determined from one complete 
orbit, on September 15 , 1996 . This orbit, which includes an episode of substantial enhancement 
of BrO in the Antarctic spring, has been selected for the BrO intercomparison being carried out 
by the various GOME researchers. 
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Figure 1. GOME BrO for overflight of EIRS-2 1997 Alaska flight 
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Figure 2a. GOME Vertical column BrO vs. solar zenith angle: Sept. 15, 1996 
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Figure 2b. GOME Vertical column BrO vs. latitude: Sept. 15, 1996 




GOME and SCIAMACHY Studies Through September 1997 


The research at the SAO for the GOME and SCIAMACHY projects has included development 
and maintenance of scientific software for the GOME retrieval algorithms, advising and as- 
sisting in the development of GOME operational software, participation in GOME validation 
and science studies, sensitivity and instrument studies to help finalize the definition of both 
instruments, completion of the SCIAMACHY study to produce the Scientific Requirements 
Document for Data and Algorithm Development, development of scientific software for GOME 
and SCIAMACHY retrievals, and consultation on the development NRT and OL operational 
data processing and data products. 


GOME 

The development of the scientific code that implements the GOME retrieval algorithms and their 
realization as the first release of operational processing code was finished. Following the end 
of the GOME commissioning phase (March 1996) and the initial validation of the operational 
retrieval algorithms, the first official release of GOME level 2 products (total ozone columns) 
took place in July 1996. GOME level 1 data products (radiance spectra) were released in late 
1997. SAO GOME studies include the following items: 


• Guiding the development of GOME level 0-1 and 1-2 operational software through exten- 
sive collaboration with the DLR. 

• Cloud studies, including investigation of fitting schemes and analysis of uncertainties due 
to various error sources, including line parameters of O 2 . These resulted in production 
of the prototype operational algorithm for GOME cloud correction which has now been 
implemented in the level 1 operational software. 

• Development and incremental improvement of the GOMEware scientific software for thor- 
ough scientific testing of GOME algorithms, including those that have been implemented 
initially in the operational framework at DLR. The activity was centered at SAO, but 
included collaboration with a number of the European institutions involved in GOME. 

• Participation in GOME validation and scientific studies. Note that the validation studies 
include ozone profiles and tropospheric ozone abundances, which are being produced as 
off-line scientific products. 



Development of fitting methods for direct (non-DOAS) fitting of GOME radiances. This 
now provides the best available global measurements of BrO from GOME (by an order 
of magnitude) and will provide substantially enhanced S 0 2 measurements (pollution and, 
possibly, background levels) in the near future. It has also permitted SAO to make the 
first and, to date, only detections of CIO from GOME measurements. This method allows 
single, non-smoothed GOME radiance spectra to be fitted to RMS levels of, typically, 
5 xl 0 “ 4 , in optical depth units. 

Collaboration on the development of the DOAS method as the operational procedure 
for GOME level 1-2 processing to produce O3 columns from GOME. This collaboration 
includes code development, measurement window selection, and a number of related issues. 

Changes to operational code for GOME level 0-2 processing, to adjust for post launch 
conditions and to correct initial deficiencies. 

Participation in GOMEtran finite-difference forward model development, in collaboration 
with the University of Bremen. 

Spectroscopic and aerosol database development, including the production of an improved 
database of molecular parameters for the visible 0 2 A band, and improved wavelength 
calibration of the BrO spectrum. 

Radiative transfer studies related to wavelength calibration and spatial and spectral alias- 
ing, including development of an improved solar reference spectrum for wavelength- specific 
applications. 

Improved determination of the wavelength-dependent Rayleigh scattering cross section. 

Improved molecular parameters for the Ring effect through molecular physics studies em- 
ploying the best currently available molecular data; improved atmospheric cross sections 
describing the Ring effect (the last three items are included in “Ring Effect Studies: 
Rayleigh Scattering, Including Molecular Parameters for Rotational Raman Scattering, 
and the Fraunhofer Spectrum,” Chance and Spurr, 1997 . 

Improved wavelength calibration of GOME irradiance and radiance spectra, through cross- 
correlation and nonlinear least-squares techniques. For example, in the BrO fitting region, 
we routinely perform wavelength calibrations of GOME data to substantially better than 
0.001 nm in absolute vacuum wavelength. 

Participation in the ESA cloud and aerosol ad hoc study group to produce cloud and 
aerosol data products from GOME. 



• Guiding and advising work on operational processor updates and reprocessing of GOME 
data. This includes, for example, the future use of the TOMS Version 7 ozone climatology 
in GOME data processing. 

• Participation in the GOME Scientific Advisory Group (GSAG) and chairing the data and 
algorithm subgroup (GDAG). 


SCIAMACHY 

SCIAMACHY is a joint German/Netherlands/Belgian program, with scientific participation 
from other European countries and, in the U.S., the SAO. It is included for launch aboard the 
ESA Envisat-1, currently planned for early 2000. 

The SCIAMACHY Science Advisory Group (SSAG) was formed, including K. Chance of the 
SAO as a U.S. member and chair of the scientific working subgroup for Algorithm Development 
and Data Usage (SADDU), J. Geary of the SAO as a U.S. member and advisor on instru- 
ment and detector issues, and R. Spurr as the third U.S. member and organizer and chair of 
the GOME/SCIAMACHY Scientific Working Sessions. SAO SCIAMACHY studies under the 
present grant include the following items: 


• Participation in the definition of SCIAMACHY instrument, including: consultation on 
SCIAMACHY detector, cooling and instrumentation issues; finalization of the SCIA- 
MACHY band definitions and detector selection criteria; synthetic retrieval studies for 
atmospheric spectroscopy in the infrared contributing to these definitions. 

• Chairing the scientific working subgroup for Algorithm Development and Data Usage; Or- 
ganizing and chairing the joint GOME/SCIAMACHY Scientific Working Sessions. There 
have been six of the to date. The seventh will be held in April 1998, and in early 1999 there 
will be an ESA-sponsored atmospheric measurements conference that is an outgrowth of 
these sessions. 

• Heading the development, compilation, and production of the Scientific Requirements 
Document for SCIAMACHY Data and Algorithm Development. 

• Contributions to the SCIAMACHY Instrument Requirements Document 

• Contributions to the SCIAMACHY Scientific Requirements Document 



• Radiative transfer studies, including development of line-by-line modeling capability for the 
infrared and visible and ray tracing for the limb and occultation measurement geometries. 
The issue of properly combining limb and nadir measurements is being addressed by the 
SADDU, in conjunction with the development of the SCIAMACHY observational strategy. 

• Development of the Basic Infrared Absorption Spectroscopy (BIAS) technique for appli- 
cation to the SCIAMACHY infrared channels. 

• Extension of GOME software development to SCIAMACHY. Existing GOME databases 
are extended in wavelength where appropriate to near infra-red. Line-by-line spectroscopic 
data sets have been compiled for SCIAMACHY (infrared) channels 7 and 8. 

• Participation in the definition of Near- Real- Time (NRT) and Off-line (OL) data prod- 
ucts for SCIAMACHY, including close collaboration with the DLR in the generation of 
SCIAMACHY NRT specification documentation. 
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ABSTRACT 

The Global Ozone Monitoring Experiment (GOME) and the SCanning Imaging Absorption spectroM- 
eter for Atmospheric CHartographY (SCIAMACHY) are diode array-based spectrometers that will make 
atmospheric constituent and aerosol measurements from European satellite platforms beginning in the mid 
1990’s. GOME measures the atmosphere in the UV and visible in nadir scanning, while SCIAMACHY 
performs a combination of nadir, limb, and occultation measurements in the UV, visible and infrared. We 
present a summary of the sensitivity studies that have been performed for SCIAMACHY measurements. 
As the GOME measurement capability is a subset of the SCIAMACHY measurement capability, the nadir, 
UV and visible portion of the studies shown here apply to GOME as well. 

1. INTRODUCTION 

SCIAMACHY sensitivity studies include its three measurement geometries: nadir observing, viewing 
the earth’s limb in scattered light, and solar and lunar occupations. In nadir observation alone, profile 
information is determinable in some cases by differential penetration of backscattered light at different 
wavelengths (as in TOMS/SBUV), and by the variable temperature structure of some molecular absorp- 
tions. Height resolution from nadir measurements is limited to 8-10 km from differential penetration and 
3 km in favorable cases using temperature structure. Height resolution in limb viewing is 3 km, limited 
by the weighting functions for limb scattering and by the spacecraft stability. Height resolution is also 
3 km for occupations, limited primarily by the telemetry data rate; 1 km resolution would otherwise be 
possible. The quantities retrieved from SCIAMACHY measurements include: 

Nadir observations: Total column amounts of O 3 , 0 4 , O 2 , CO, H 2 O, CH 4 , CO 2 , NO 2 , N 2 O, HCHO, 
S0 2 , CIO, OCIO, and BrO; stratospheric profiles of 0 3> CO, H 2 0, CH 4 , C0 2 , N 2 0; the column of NO 
above the ozone layer; tropospheric profiles or columns of O 3 , CO, H 2 O, CH 4 , CO 2 , N 2 O. Stratospheric 
profile information, discrimination between stratospheric and tropospheric columns, and, in some cases, 
tropospheric profile information is derived from the temperature dependences of the absorption features. 
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Limb viewing observations; Stratospheric profiles of O* (20-50 km), 0 2 (*A) (50 km-90 km), 0 2 
(20-50+ km), CO ( 20-35 km), H 2 0 (20-53 km), CH 4 (20-40 km), C0 2 (20-50+ km), NO (40-140 km), 
NO* (20-40 km), N 2 0 (20-35 km), BrO (20-25 km under normal conditions; 20-30 km under ozone hole 
conditions). 

Solar /lunar occult at ions; Stratospheric profiles of O 3 (15-50 km), 0 2 (15-50+ km), CO (15-35 km), 
H 2 0 (15-53 km), CH 4 (15-40+ km), C0 2 (15-50+ km), N0 2 (15-40 km), N 2 0 (15-35 km), BrO (15-35 km 
under normal concentrations). Occult&tions are more infrequent than limb scattering observations, and 
are limited to latitudes of 50-80°, but provide the most sensitive and precise profile retrievals. 

Tropospheric amounts from limb-nadir subtraction: TVopospheric columns of O 3 , CO, H 2 0, CH 4| 

co 2 , N 0 2 , n 2 o. 

Clouds and aerosols: SCIAMACHY will obtain aerosol data in occultation, limb-scattering and nadir 
measurements. Stratospheric aerosol measurements include background aerosol, Junge layer, and polar 
stratospheric clouds (PSCs); tropospheric aerosols, including clouds, sand winds, and soot from forest fires. 

2. MOLECULE SENSITIVITY STUDIES 

The sensitivities for retrievals of molecular species concentrations as functions of altitude have been 
extensively studied, under a wide variety of geographic conditions, in order to accurately gauge the capa- 
bilities of the SCIAMACHY instrument. Calculations are performed for all of the SCIAMACHY modes of 
measurement, including nadir viewing, viewing of the earth’s limb in scattered light, and solar and lunar 
occultations. Note that simultaneous nadir and limb measurements are the normal mode of operation for 
the majority of dayside portion of each orbit. The calculations include realistic expectations of detector 
quantum efficiency and noise performance, and the instrument etendue and optical throughput that are 
derived from the Phase A instrument study. 

Retrieval studies for atmospheric molecules are performed using atmospheric concentration profiles from 
the MPI Mainz global 1-D and 2-D models (C. Bruhl, E. Hennig and P. Crutzen, private communication), 
calculated with clean and polluted air, daytime, conditions. Additional information on nighttime NO 3 
comes from Norton and Noxon, 1986. 1 Constant mixing ratios are assumed for C0 2 (340 ppm) and 0 2 
(20.95%). The simulations assume exposure times of 2 seconds in nadir spectra, 1.5 second in limb 
spectra and lunar occultations (appropriate to 3 km vertical resolution), and short exposure times in 
solar occultation, with the constraint that exposures are limited to 10* photons/pixel in the most highly 
illuminated portion of each detector array. Interferences from other absorbing species are included in the 
studies for each individual molecule. Current studies are limited to a maximum pressure altitude of 60 km 
because of the limitation of the model calculations to the stratosphere. In practice, profiles for at least 
0 2 , O 3 , and C0 2 will extend to higher altitudes. 

Spectra for the simulations are taken from the best available sources, including digitally acquired electronic 
spectra obtained in the laboratory at MPI Mainz. Spectra in the infrared portion are simulated using the 
1986 HITRAN line parameter listing . 3 

Calculations of geometries are performed with multilayer, curved-shell models, including refraction, to 
provide line of sight column densities as functions of temperature, pressure, and altitude. Calculations in 
the UV and visible parts of the spectrum are made using the AFGL LOWTRAN7 code, including multiple 
scattering in the nadir and limb geometries. Calculations in the infrared are performed with LOWTRAN7 
and with the Smithsonian Astrophysics! Observatory line-by-line radiance code . 3 Most calculations are 
made with the earth’s average albedo of 0.3. The effect of surface spectral reflectance on nadir measure- 
ments was investigated in more detail, as discussed further below. Lunar occultations are considered for 
the measurement of NO 3 , as discussed in the section for that molecule. 
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The 2-D model results used in the calculations are for January 1, with global coverage extending from 85° S 
to 85° N, which samples a very representative set of the earth’s atmospheric conditions. Most calculations 
are done at 35°N, which corresponds to a solar zenith angle of 63.4° for the EPOP 10:00 crossing time. 
The effect of other solar zenith angles on measurements in the nadir and limb geometries, including 
measurements in polar regions, was investigated in some detail, as discussed below. 

The solar illumination used in these studies is the LOWTRAN7 solar source function. Aerosol and molec- 
ular scattering are also from LOWTRAN7. 20 km is taken as a practical lower limit for limb scattering 
measurements due to confusion by multiple scattering at lower tangent heights. 15 km is taken as the 
lower limit for occultation measurements because of atmospheric refraction and extinction by aerosols. It is 
hoped that, in practice, both of these limits can be extended downward. The anticipated EPOP spacecraft 
jitter of 0.01° corresponds to 0.57 km at the limb, or 8% in airmass. We combine this with the effect of 
multiple scattering on limb retrievals: Mount et uL, measured NO 2 in limb scattering with negligible 
effects due to multiple scattering down to 28 km, with the Solar Mesosphere Explorer satellite, using wave- 
lengths from 430 to 450 nm. 4 SC1AMACHY obtains significant information on NO 2 and O 3 out to 500 nm 
implying, for Rayleigh scattering varying as 1/A 4 , negligible confusion due to multiple scattering down to 
25 km. Thus, we estimate limb measurement errors in precision for profile determinations, discussed fur- 
ther below for individual molecules, of 10% above 25 km and 20% from 20-25 km, unless further limited by 
the signal to noise ratios of the measurements for such molecules. Determination of stratospheric columns 
from limb measurements in this wavelength range, particularly for species with concentrations which peak 
in the stratosphere, will be somewhat more precise, since the multiple scattering predominantly causes 
smearing of profile information over a scale height about the tangent altitude. Species measurements in 
the stratosphere using limb scattering at significantly shorter wavelengths will provide only limited profile 
information due to the larger effects of multiple scattering. Limb measurements at longer wavelengths, 
particularly in the infrared, are less affected by multiple scattering. Nadir and solar occultation measure- 
ment capabilities are calculated using a lower limit of 1% for measurement precision due to the difficulty 
of making more precise measurements and calibrations. We hope, of course, to be able to achieve higher 
precision in practice so that this technique will be capable of determining long-term trends in global 0 3 - 

Sensitivity and interference studies for some species to be measured in the infrared are checked by examining 
a high resolution (0.007 cm* 1 ) solar spectrum from Kitt Peak National Observatory (KPNO), obtained 
courtesy of D. W. Johnson and G. Stokes. The spectrum was taken as part of the BIG campaign, under 
clean tropospheric conditions, at 1.1 airmass. Confirmation of O 4 is from a 1.0 airmass high resolution 
visible solar spectrum from KPNO (courtesy of R. Kurucz). 

Figure 1 shows the solar spectrum as measured by SCLAM ACH Y , including the etendue, optical through- 
put, detector quantum efficiency (normalized to photons pixel 1 s 1 for the various detector chan n els), 
and the beam dilution (18% of the instrument field of view is filled by the sun when SCLAM ACHY views 
the central part of the solar disk) . 

A molecule by molecule discussion of the retrieval studies follows. The results are summarized in Figure 
2, which shows altitude ranges for profile retrievals, and Table 1. The altitude limits shown here are for 
retrieval to a 3 km vertical resolution with <10% uncertainty in the concentrations, excluding uncertainties 
in cross-sections, for the stratosphere, plus the altitudes over which any tropospheric information can be 
obtained. These limits are derived from the January 1, 35°N studies, implying a solar zenith angle of 
63.4°, as mentioned previously. For most molecules they are valid for solar zenith angles of up to at least 
80°. Molecules which may possibly be measured by SCLAM ACHY but are not included in the the studies 
due to the lack of suitable available laboratory spectra include 0( 3 P), HO 2 * and HONO. 

2.1 Molecules in the normal atmosphere 

O3 concentrations are retrieved from SCIAMACHY measurements using the ultraviolet (Hartley and 
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FIGURE 2. SCIAMACHY stratospheric and tropospheric profile ranges 
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Table 1 


Quantities Retrieved from SCIAMACHY Nadir /Limb Observations 


Species 

Retrievable Quantity* 

Wavelength (nm) 

Notes/ Applications 

o 3 

Profile (S, T) 

255-350, 480-680 


Oj 

Profile (S, T) 

690, 760 

Cloud tops/ 
boundary layer 

Oj( l A) 

Profile (S) 

1260-1275 

Mesospheric O 3 

0 4 

Column (T) 

475, 530, 560, 630 

Cloud tops/ 
boundary layer 

h 2 o 

Profile (S, T) 

700-1000, 2250-2390 


CO 

Profile (S), Column (T) 

2300-2390 

Limited tropospheric 
profile information 

COa 

Profile (S, T) 

1980-2020 

P, T retrievals 

ch 4 

Profile (S, T) 

2250-2390 


HCHO 

Column (T)t 

310-360 

Biomass burning, 
tropical vegetation 

S0 2 

Column (T)t 

290-310 

Industrial pollution, 
volcanos 

NO 

Profile (S) 

255-280 

Above 40 km 
in emission 

NOj 

Profile (S), Column (T)t 

300-600 

Lightning, 

combustion 

NO 3 

Profile (S) 

660-665 

Lunar occultation 

n 2 o 

Profile (S), Column (T) 

2250-2300 

Limited tropospheric 
profile information 

CIO 

Column (S)^ 

300-310 

Polar regions in 
spring 

0C10 

Column (S)^ 

320-420 

Polar regions in 
spring 

BrO 

Profile (S) 

310-345 



• S = stratosphere; T = troposphere. Mesosphere not included in this table, 
t Observable in regions with relatively high concentrations, 
t Observable in perturbed “03 hole” regions. 
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Huggins) bands and the visible (Chappuis) bands. Precise determinations of column O3 are obtainable 
from Chappuis band measurements. Measurements of the UV bands, using knowledge of penetration 
depth as a function of wavelength to obtain stratospheric profile information, has become a standard 
technique for satellite ozone measurements. This differential penetration is measured using a number of 
discrete wavelength channels that are much broader than the SCIAMACHY resolution elements. Thus, 
SCIAMACHY includes this type of measurement as a subset of its O3 measurements. 

SCIAMACHY can also measure ozone profiles in the limb scattering geometry, up to at least 50 km 
in most geographic locations. Although the signal to noise ratio in these measurements is quite high, the 
precision of concentration profiles from limb measurements will be limited to ca. 10% above 25 km and 
20% from 20-25 km due to the mixing of light from various sources in the multiple scattering process. 
Stratospheric column O3 from limb measurements can be determined to 10% since limb measurements 
extend down to at least 1/2 scale height below the O3 peak. 

A strength of the SCIAMACHY instrument is that it can provide extensive maps of tropospheric 
ozone. Some determinations of tropospheric O3 in tropical regions have been made in recent studies by 
subtracting stratospheric O3, measured by the SAGE instrument, from TOMS total ozone measurements. 6 
SCIAMACHY includes the UV bands used in past and present satellite instruments for measuring O3, 
but affords continuous coverage of the region above 240 nm, and at significantly higher resolution. This 
allows us tremendous leverage in retrieval of tropospheric ozone profiles due to the temperature dependent 
structure of the Huggins bands. 6,7 The Huggins bands have discrete vibrational structure between 300 
and 370 nm, with features having widths significantly less than 1 nm. This is the sharpest vibrational 
structure of any electronic band of O3 that has been studied. This structure has strong temperature 
dependence due to the onset of thermal population of excited vibrational levels in the electronic ground 
state. It is this temperature dependent structure that provides a spectral signature for tropospheric ozone; 
its differential character derives from the onset of weaker band absorption between the stronger absorption 
peaks, which are due to absorption from O3 (000). One may think of the differential of the O3 spectrum 
with temperature as part of a basis set used to match the observed spectrum. Note that tropospheric 
ozone is significantly warmer than stratospheric ozone at any altitude in the stratosphere where there 
is significant O3. Thus, the tropospheric O3 column is simply obtainable from an instrument with the 
resolution and sensitivity of SCIAMACHY. Determination of the complete altitude profile, including the 
troposphere, requires the development of a more sophisticated algorithm, which will include differential 
penetration and the Chappuis band column measurements. Figure 3 shows the sensitivity of the nadir 
spectrum as measured by SCIAMACHY, at the SCIAMACHY resolution, for 1% changes in total ozone 
column at various levels in the atmosphere. The sensitivity is shown as signal to noise ratio per pixel 
that can be achieved in 2 seconds of averaging time by SCIAMACHY in measuring these 1% differences at 
different altitudes. Figure 3 was performed using the full multiple scattering formalism of LOWTRAN7, to 
show the larger characteristics of the radiance changes due to ozone distribution. This does not, however, 
fully model the temperature dependence of the bands. Our sensitivity studies indicate that SCIAMACHY 
can achieve a precision corresponding to at least 1% of the total column ozone for retrieval of O3 to 3 km 
vertical resolution throughout the stratosphere and troposphere. 

We have calculated the ability of SCIAMACHY to measure polar ozone depletion. Calculations are 
performed using a solar zenith angle of 80°, appropriate to 80° N-S at the equinox. Moderate ozone 
hole conditions are simulated by removing 20% of the total column ozone from atmospheric above the 
tropopause, symmetric about the ozone peak. The nadir measurements will be able to measure polar 
ozone depletion to the full geographic coverage of the orbit. SCIAMACHY should also be able measure 
polar ozone depletion in the limb geometry to even higher latitudes, down to quite low solar zenith angles. 
The geographic limits to high latitude for this type of measurement are difficult to estimate precisely at 
present due to the difficulty in estimating the light levels under these conditions. 

As an additional check on the ability of SCIAMACHY to retrieve stratospheric ozone profiles, and to 
distinguish between stratospheric and tropospheric ozone, we have calculated the ability of SCIAMACHY 
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Figure 3. SC1AMACHY S/N for measurement of 1% 0 3 changes (nm) 




to measure stratospheric O3 in solar occultation, chiefly using the Chappuis bands, so that the tropospheric 
residual may be obtained by subtraction. The signal to noise ratio per pixel for measurements of the 
Chappuis bands at their peak for these measurements, which correspond to a vertical resolution limited 
only by the slit size (1 km), is tremendous, and there are a large number of pixels covering the Chappuis 
bands. The implied measurement precisions are extremely high, and will be limited in practice by other 
factors, such as spacecraft pointing and intensity calibration. SCIAMACHY will be able to measure O3 
to better than 1% precision for each 3 km height interval from 20-40 km, assuming solar tracking to 
sufficient precision, and to better than 5% for the total stratospheric column. This is sufficient precision to 
quantitatively test the determination of tropospheric O3 from the temperature dependence of the Huggins 
bands under conditions of relatively high tropospheric 0 3 , and to provide for some independent mapping 
of tropospheric O3 at modest precision and geographic coverage. Stratospheric measurements in solar 
occultation extend at lower precision to at least the stratopause. 

In summary, SCIAMACHY can determine the ozone column density quite accurately from the nadir 
measurements (to better than 1% precision); it can determine profiles to better than SBUV precision 
from the nadir measurements; it can cleanly distinguish between tropospheric O3 and stratospheric O3 
from the nadir measurements and give moderate (1% of total column ozone) profile information in the 
troposphere; it can provide further stratospheric profile information at moderate precision from the limb 
scattered light; it can provide stratospheric profile information (and possibly upper tropospheric profile 
information) with very limited geographic coverage from solar occultation measurements, with sufficient 
precision to test the algorithm for determination of tropospheric O3 from nadir measurements; it can 
make precise measurements of polar ozone depletion to at least 80° latitude, and probably further. The 
SCIAMACHY measurements of ozone are not subject to long-term degradation of the type suffered by 
SBUV/TOMS, as discussed in the instrumentation section of this report. 

O a ( 1 A). Quantitative laboratory Bpectra of Oaf 1 A) are not available for a sensitivity limit study, but the 
positions are well known. Absorption by tropospheric ground state oxygen to the 1 A) state is quite promi- 
nent in the infrared spectrum. This absorption masks the emission by Oaf 1 A) in the nadir observations. 
Oa( x A) is used by the Solar Mesosphere Explorer satellite in limb scattering observations to determine 
the O3 profile from 50-90 km; 8 SCIAMACHY will be able to measure it in limb scattering for the same 
purpose. 

O4. Quantitative laboratory 0< spectra are not available for a sensitivity study. As it is quite prominent 
in the solar spectrum observed from the ground. As it is a van der Waals molecule, and thus has a 
concentration strongly dependent upon temperature and pressure, it is a useful diagnostic for cloud top 
and boundary layer height. 

O2 concentrations provide a cross-check of the pointing calibration in occultation and limb measurements. 
In the nadir view, penetration depth determinations from 0 2 provide information about cloud and bound- 
ary layers heights. O2 absorbs very strongly in the atmosphere, chiefly near 690 and 760 nm. 0 2 column 
measurements will have a precision of 1%, our projected limit for strong absorbers. 0 2 measurements in 
the limb and occultation geometries will also have the limiting sensitivity up to at least the stratopause. 

CO. The 2 < 0 vibrational overtone band of CO lies between 2.3 and 2.4 fan. SCIAMACHY can be used in 

nadir view to measure the CO column abundance with a sensitivity limit of 2 x 10 16 cm 2 , corresponding 
to 1% of background tropospheric CO, when the earth’s reflectance at these wavelengths is as low as 
0.1. The simulations have taken proper account of solar CO, which has lines of comparable size to the 
absorption from background terrestrial CO at about m = 9 and higher, but which is divided out by our solar 
calibrations and which, in any case, has very characteristic spectrum differences due to the much higher 
temperatures of solar CO. Since the tropospheric CO measurements are fundamental to the success of 
SCIAMACHY, the results of the simulation are checked by examining the KPNO infrared solar spectrum. 
The CO lines are quite prominent in the spectrum, as predicted, with most lines being sufficiently separated 
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from interfering HjO and CH 4 to be cleanly measurable at the SCIAMACHY resolution. Additionally, 
** bave calculated the ability of SOIAMACHY to measure CO in the limb geometry, by comparing the 
limb-scattered radiation strength with the MPI model CO concentration values and the CO line strengths. 
Stratospheric CO absorption can be measured with a precision ranging from 1.5% at 20 km to 15% at 
30 km (by 35 km the precision has dropped to 50%). In practice, due to effects of multiple scattering, 
the determination of total stratospheric CO will be precise to 10-15%. Since only ca. 5% of the column 
CO is expected to be in the stratosphere, the limb measurements will provide for a correction to the 
column tropospheric CO to better than 1%. Limited CO profile information for the troposphere will be 
available from the differential temperature dependence of the various CO lines in nadir observing. CO is 
also measurable at very high signal to noise ratios in solar occultation spectra. 13 CO has line strengths 
down by a factor of 0.011 from **CO, due to the 13 C isotopic abundance. It will be readily measurable in 
the lower stratosphere in solar occultation. There is a marginal possibility that the column of 13 CO will 
be measurable from signal averaging of nadir spectra. 

HjO absorbs strongly between 700 and 1000 nm, and in the vibrational overtone bands near 2400 nm. The 
strong temperature dependence of these latter bands allows retrieval of the altitude profile of tropospheric 
HjO in the nadir view. HjO is retrievable in all SCIAMACHY measurement geometries to the limits 
of the SCIAMACHY technique - 1% in nadir, 10-20% in stratospheric profile from limb measurements, 
1% in stratospheric profile from solar occultation, and 10% per 3 km interval in tropospheric profile from 
deconvolution of the temperature-dependent nadir spectrum. 

CH 4 has strong vibration overtone absorption from 2200 to 2400 nm. As with H 2 0, the temperature 
dependence allows retrieval of the tropospheric altitude profile from the nadir view. The limits for CH 4 
profile measurements in limb-scattered light extend to 40 km in the stratosphere, due to the high concen- 
tration of CH 4 and the strength of the absorption. The column amount of 13 CH 4 will almost certainly be 
measurable with SCIAMACHY to good precision, since 1.1% of 12 CH 4 is readily measurable. However, 
line positions and strengths for 13 CH 4 are not currently catalogued, so it is not possible to further quantify 
its measurement. 

COj has its strongest overtone absorption between 1900 and 2200 nm. As with H 2 0 and CH 4 , the 
temperature dependence can be used to retrieve the tropospheric profile in the nadir view. The C0 2 
absorption is so strong that the retrieval of C0 2 in all geometries is limited only by the SCIAMACHY 
baseline precision considerations. We note here that the SCIAMACHY C0 2 measurements are fully as 
capable of providing atmospheric pressure and temperature information as the standard 15 pm infr ared 
technique, since the band used (the strongest overtone absorption of C0 2 ) originates from the ground 
state and thus supplies the same information about the distribution in energy levels of atmospheric C0 2 . 
It may be possible to use this C0 2 information to partly correct for the effects of spacecraft jitter on limh 
measurements, but such corrections have not been included in the present calculations. 

NO. The NO profile above the ozone layer (from 40-140 km) can be determined from SCIAMACHY limb 
and occultation measurements using the (1,4), (1,6), and (0,3) gamma bands which fluoresce from 255-280 
nm, as determined in model sensitivity studies by Frederick and Abrams. 9 In nadir view, gamma band 
emission gives a column amount for NO above where O3 absorbs strongly. 

N0 2 concentrations are retrieved from the uv- visible absorption between 300 and 600 nm. In limb viewing 
NOj is readily observed in the stratosphere as already shown by SME, and by balloon measurements using 
this technique. 10 Stratospheric N0 2 is measurable to 10% precision from 25-40 km and 20% precision from 
20-25 km in limb viewing, limited by the considerations discussed earlier. Stratospheric column determi- 
nations from limb viewing are good to 10%. In occultation, the profile will be measurable to 1% precision 
from 40 to below 20 km. In nadir viewing the total column is observed to the limiting 1% precision. From 
these two measurements tropospheric N0 2 amounts may be determined under some conditions. SME 
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measured N0 2 column amounts above 28 km in the range 1-3 xlO 15 molecule cm - ’. 4 ' 11 This corresponds 
closely to the total N0 2 column abundances measured by SAGE . 12 A background tropospheric mixing 
ratio of 30 pptv yields a tropospheric column amount which is approximately an order of magnitude lower. 
In background clean sir conditions deconvoluting the stratospheric and tropospheric amounts from SCIA- 
MACHY data is marginal. In polluted air this task becomes much simpler. For 3x 10 15 cm " 2 stratospheric 
N0 2 , the column uncertainty is 3xlO u cm -2 , which corresponds to a tropospheric limiting concentration 
of 16 pptv for S/N = 1. Tropospheric N0 2 in even moderately polluted regions is substantially larger 
than this. In such regions, the limi t for tropospheric column becomes the 1% precision of the total column 
N0 2 . The tropospheric column can be determined to at least 10% precision for concentrations above 160 
pptv. The N0 2 absorption features are known to be temperature dependent, which in nadir view may 
assist further in the differentiation of tropospheric and stratospheric amounts. 

NO 3 . The nighttime concentration of NO 3 has been measured successfully by absorption of lunar and 
stellar light . 13 The strongest N0 3 absorption, between 660 and 665 nm, is 2-3% deep in the occultation 
geometry for tangent heights of 28-38 km. Using lunar occultation we will be able to measure the profile 
of NO 3 , when the lunar cycle permits, with potential long-term latitude coverage of 50-80° N and S . 14 
The scan system can observe in both hemispheres, but the spacecraft accommodation will probably imply 
that the moon will only be observable in one hemisphere. Studies will be made in Phase B of the possible 
use of planets and stars (e.g., Venus) for nighttime occultation measurements. This is a significant portion 
of the region where enhanced NO 3 is expected in polar night conditions. The lunar occultation signal 
measured by SCIAMACHY is down from the solar occultation signal by a factor of 1.2 xlO " 6 due to the 
geometric factors and the lunar albedo. For N0 3 at 30 km, using 1.5 second integration time, the numbers 
of photons per pixel at 660-665 nm is 4.2 x 10 6 for a signal to noise ratio per pixel of 2.0 x 10 3 . For 2% 
absorption at band center, the NO 3 is measured with 2.4% precision. 

N 2 0. The 00°2 ♦- 00°0 vibrational overtone band of N 2 0 lies between 2.25 and 2.30 /im. Lines of N 2 0 are 
quite prominent in the Kitt Peak solar spectra, with many individual lines well separated from interferences 
at the SCIAMACHY level. The N 2 0 column is measurable to the limiting nadir precision (1%) in all but 
the lowest surface reflectance and highest solar zenith angle conditions. In limb measurement, N 2 0 can 
be measured to a precision limited only by spacecraft jitter (< 10%) up to 30 km, and to about 40% up 
to 35 km. The total stratospheric N 2 0 will be precise to 10%; tropospheric columns should be obtainable 
by subtraction to 5% or better. The temperature dependence of the multiple N 2 0 lines may provide some 
information on the vertical distribution of tropospheric N 2 0 . 

2.2 Polluted troposphere molecules 

HCHO absorbs between 250 and 360 nm, with its strongest absorption between 300 and 350. The short 
wavelength absorption is strongly obscured by O3 absorption, leaving a contribution between about 300 
and 350 nm. This window is highly dependent on the 63 concentration and the viewing geometry and 
the temperature and concentration of tropospheric O3 (increased absorption by warm tropospheric ozone 
will further obscure the HCHO). HCHO is measurable in the nadir view, chiefly at low latitudes. The 
sensitivity for HCHO measurements with this concentration is 5. Thus, moderately enhanced HCHO can 
be observed with modest precision and major pollution can be readily measured. Note that this signal to 
noise ratio is determined using only the portions of the HCHO spectrum close to 350 nm, which remain 
reasonably unobscured by 0 3 . The sensitivity increases by as much as a factor of 2 under conditions of 
low column O3. 

S0 2 absorbs between 290 and 315 nm. It is observable in the nadir view, and under conditions of low 
tropospheric ozone, in the polluted troposphere. For a tropospheric concentration of 20 ppbv (moderate 
pollution in the boundary layer) SCIAMACHY achieves a signal to noise ratio of 60 (i.e., better than 2% 
p recision ), but this is highly variable, depending on the O3 concentration. S0 2 in normal background 
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concentrations is not observable. 

2.3 Ozone hole molecules 

CIO absorbs from 220 to 310 nm. It is strongly masked by ozone absorption. Our sensitivity limit for 
the nadir view varies from about 10 13 to 10 15 cm -3 , depending on the ozone concentration. We do not 
expect to be able to detect CIO under normal conditions. Unusually high CIO concentration coupled with 
low O3 concentration (which will improve the detection limit substantially), as found in the Antarctic 
ozone “hole” 15 will make the CIO column observable. An absorption spectrum of CIO calculated using a 
representative 0 3 hole concentration of 1.5xl0 15 cm -3 , compared with a synthesis of the nadir spectrum 
at an 80° solar zenith angle, under moderate ozone hole conditions, gives a precision for CIO column 
measurements of 10%. 

OCIO absorbs strongly between 280 and 440 nm. It is formed in the atmosphere from the reaction of BrO 
with CIO and destroyed by photolysis. OCIO has been observed using its structured absorption between 
320 and 420 nm at night over Antarctica in spring (column densities of lxlO 14 molecule cm -3 . 16 - 17 Using 
our synthesized nadir spectrum for moderate ozone hole condition and the OCIO absorption for lxlO 14 
cm -3 , corresponding to a twilight concentration, we obtain a 2% measurement precision for OCIO. For a 
typical daytime concentration (2 xl0 ls cm -3 ) we can measure to 10% precision. It may be possible to 
obtain OCIO profile information at ozone hole concentrations using lunar occultation; insufficient profile 
modeling is available to accurately estimate the profile limits. 

BrO has been measured in the Antarctic spring using its UV absorption. 17 By comparing the absorption 
of BrO for a nadir concentration of 2.5 xlO 13 for our ozone hole geometry, we determine that we can 
measure this column density of BrO with a precision of 10%. Under normal conditions, the BrO column 
density is 1.5 xlO 13 , implying measurements to 15%. BrO profile ranges are estimated using MPI 2-D 
models calculations. With solar occultation, BrO in normal concentrations can be measured to 25% from 
15-35 km. In limb measurements, BrO profiles can be determined to 50% precision from 20-25 km in 
normal concentrations and ca. 30% in ozone hole concentrations. 

2.4 Surface reflectance and solar zenith angle studies 

Absolute radiometric calibration of the SCIAMACHY instrument, combined with direct solar observations, 
will provide new global information about surface spectral reflectance in the range 300-2400 nm. Although 
the spatial resolution of SCIAMACHY is limited, this data will be very useful scientifically. Further studies 
of the retrieval of the reflectance measurements by SCIAMACHY will be made in Phase B. 

Surface reflectance and solar zenith angle considerations for molecule measurements in nadir viewing are 
coupled together since the increased path due to larger solar zenith angle can offset decreased light due 
to lower reflectance for those molecules measured at wavelengths that penetrate fully to the ground. The 
increased path can also exacerbate measurements of tropospheric molecules made at shorter wavelengths 
due to increased loss of tropospheric information because of increased Rayleigh scattering in the strato- 
sphere. Our standard conditions for sensitivity studies (35° N: 63° solar zenith angle; albedo =0.3) are 
fairly stringent in the sense that, for the suite of molecular measurements by SCIAMACHY, the deterio- 
ration in global coverage due to lower reflectance and other solar zenith angles is not severe. The losses 
that do occur are mainly in tropospheric measurements. 

Surface reflectances and solar zenith angles other than our standard conditions have their greatest impacts 
on the following measurements: Tropospheric Os- At 320 nm, 13% of the light measured by SCIAMACHY 
at our standard 35° N (63* solar zenith angle) is reflected from the ground. For albedo = 0.1, this decreases 
to 3.6%, and the precision for tropospheric 0 3 is a factor of 2 lower. Losses in tropospheric O3 information 
are partly compensated for by increased solar zenith angle, since the light at the red end of the Huggins 
bands has significant penetration to quite high solar zenith angle. Tropospheric NO2 is measured at 
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modest precision in our standard conditions, in a wavelength range where Rayleigh scattering is not 
a major problem. Sensitivity to tropospheric NO 2 varies roughly as the square root of the reflectance. 
Measurements of the tropospheric columns of HCHO and SO 2 will suffer with decreased reflectance. HCHO 
sensitivity will increase almost linearly with the increased path due to higher zenith angle, up to quite 
high angle, because of the long wavelength extent of the spectrum. SO 2 sensitivity changes with zenith 
angle will be small. CIO, OCIO, and BrO column measurements were calculated at a high solar zenith 
angle (80°), appropriate to measurements by SCIAMACHY in ozone hole conditions. CIO is measured at 
such short wavelengths that changes in surface reflectance are not a serious consideration in measurement 
sensitivity; variability in O 3 has a much more serious impact. CIO will generally be seen at high solar zenith 
angles, so that sensitivity vs. zenith angle is not an important issue. OCIO is also generally seen at these 
high zenith angles. Its measurement is more sensitive to surface reflectance since the spectrum extends 
well beyond 400 nm. Fortunately, the climatology of OCIO is such that measurements will generally occur 
in regions of high surface reflectance; exceptions occur when the ozone hole conditions extend over ice-free 
ocean waters. Even in these conditions, the strength of the absorption in the 300-350 nm range ensures 
useful measurement sensitivity. The spectrum of BrO has its major absorption in the 300-350 nm range, 
long enough that O 3 absorption is not a serious problem and short enough that lower surface reflectance 
also will not seriously hamper its measurement. BrO measurement under normal conditions (i.e., away 
from 03 -depleted polar regions) will be decrease as the path length decreases due to lower solar zenith 
angle. 


3. AEROSOL STUDIES 


The presence of aerosols in climatic models leads to a more important energy deposition in the stratosphere 
and to a cooling of the troposphere. More recently, the Antarctic ozone hole phenomenon has been 
explained by the conversion of H Cl to CIO by heterogeneous reactions on the surfaces of polar stratospheric 
clouds. 

The SCIAMACHY instrument will obtain aerosol data in occultation, limb-scattering and nadir mea- 
surements. The high resolution permits distinction between gases and aerosols and will allow altitude 
distributions and granulometry of the middle-atmospheric aerosols to be obtained. Proper treatment of 
aerosols in the inversion process also results in more accurate data on the vertical distribution of gases. 
Simulation and verification of the various types of spectral signatures from different aerosols are an integral 
part of the SCIAMACHY data system development plan. 

3.1 Stratospheric Aerosols 

Stratospheric aerosols have been extensively studied from balloon and space experiments over the last 30 
years: The measurements show spherical particles in the stratosphere composed of a mixture of sulfuric 
acid and water. These are divided between layers of large particle size, the Junge layer (typical size 0.15 
/im), and a background layer of aerosols of significantly smaller size (ca. 0.05 p m). However, in the polar 
stratospheric clouds associated with the Antarctic ozone hole, particle sizes larger than 1 pm have been 
conjectured. Other types of suspended particles exist in the upper atmosphere, from meteoric decay; little 
is known about them except that their metallic cores give them an important imaginary refractive index 
and that the resulting increase of Mie scattering efficiencies makes them more observable, as published in 
visual space observations and balloon data reports. SCIAMACHY, through measurements of these latter 
particles, may help solve the present controversy concerning the geophysical impacts of small meteorites 
and comets. 
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3.2 Tropospheric Aerosols 

The high spectral resolution of the SCIAMACHY instrument combined with its spatial resolution can lead 
to study of image contrast, and to a characterization of the nature of tropospheric haze in terms of natural 
phenomena. For example, sand winds have never been monitored on a global scale. Since they occur by 
definition mostly in dry regions, the absence of clouds will make SCIAMACHY a perfect sensor for them. 
SCIAMACHY will also measure the high quantities of soot from forest fires to which, for decades, the blue 
moon phenomenon has been attributed. Identification of these “clouds” and their continuous observation 
will lead to an a ss ess ment of the effects of these largely local phenomena on the global atmospheric system. 
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ABSTRACT 

The potential of the ESA Global Ozone Monitor- 
ing Experiment (GOME) to produce ozone profile 
information, has been examined by carrying out two 
sample retrievals using simulated GOME data. The 
first retrieval examines the potential of the GOME 
instrument to produce stratospheric ozone profiles 
using the traditional back-scatter ultraviolet tech- 
nique, while the second examines the possibility of 
obtaining tropospheric profile information, and im- 
proving the quality of the stratospheric profile re- 
trievals, by exploiting the temperature dependence 
of the ozone Huggins bands. 

INTRODUCTION 

The GOME instrument is a nadir- viewing grating 
spectrometer, under development by the European 
Space Agency, for launch on the second European 
Research Satellite (ERS-2) in late 1994. 

It will measure solar radiation back-scattered 
from the earth’s atmosphere in the spectral region 
240-790 nm, at moderate resolution (0.2-0.4nm) 
throughout the range. The primary measure- 
ment technique to be exploited by the GOME 
is known as Differential Optical Absorption Spec- 
troscopy (DOAS). Using this technique, the col- 
umn abundance of a trace gas is derived from a 
back-scattered ultraviolet /visible spectrum by ratio- 
ing it to a solar spectrum, high pass filtering, cross- 
correlating with a laboratory measured spectrum to 


obtain the slant column density, and converting from 
slant to vertical columns by use of air-mass factors 
(Platt & Perner, 1985; Platt et al , 1980a, 1980b, 
1981; Hubler et a/., 1984). 

In this capacity the GOME is expected to produce 
column abundances of a number of trace gases (e.g. 
0 3 , NO, N0 2 , BrO, H 2 0, 0 2 and 0 4 ). 

In addition, the GOME will produce ozone pro- 
files using the back-scatter ultraviolet method which 
has been successfully exploited by NASA’s Solar 
Backscatter UltraViolet (SBUV) instruments (Heath 
et a/., 1975). This technique requires an absolute ra- 
diometric calibration. For this purpose, the GOME 
will have the capability to observe the sun via a dif- 
fuser plate, and the moon directly. 

As the GOME will have a much broader spec- 
tral coverage and higher spectral resolution than the 
SBUV instruments, which obtain useful measure- 
ments at only eleven wavelengths between 273 and 
340 nm, it is hoped that improved vertical profile in- 
formation will be obtained. Of particular interest is 
the possibility of retrieving tropospheric profile infor- 
mation. In order to obtain ozone profile information 
in the troposphere from the GOME measurements, 
it is necessary to take advantage of the temperature 
dependence of the ozone absorption features from 
300-340 nm, known as the Huggins bands (fig. 1). 
Since the variation in absorption cross-section with 
temperature is small (around a few percent per 10K; 
Paur & Bass, 1985; Yoshino, 1988), it is necessary 
to obtain measurements with a high signal to noise 
ratio to allow the height information to be extracted. 
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Figure 1. Tie ozone absorption bands known as the 
Huggins bauds (Nicolet, 1980). 


It is expected that the GOME will be able to pro- 
vide measurements with sufficiently high signal to 
noise ratios (—1000) in this spectral region to allow 
tropospheric profile information to be obtained. The 
remainder of this paper is devoted to a discussion of 
ozone profile retrievals carried out using simulated 
GOME data. 

METHOD 

Two sample ozone profile retrievals have been car- 
ried out using simulated GOME data at a spectral 
resolution of 0.2 nm; the first using measurements 
from the 255-300 nm region only, and the second us- 
ing measurements from the 255— 340 nm region in- 
cluding the temperature dependent Huggins bands. 

The measurements were simulated using LOW- 
TRAN7, a radiative transfer code including multi- 
ple scattering (Kneizys et a/., 1988). Representative 
GOME instrument noise was then added. The simu- 
lated back-scattered spectrum was subsequently di- 
vided by the solar source spectrum in order to remove 
the effect of the large variation in incident intensity. 
The solar spectrum was also calculated using LOW- 
TRAN7. 

To date, no attempt has been made to examine the 
effect of errors in wavelength registration between the 
measured back-scattered and solar spectrum. Nei- 
ther has the effect of uncertainties in radiometric 
calibration been examined. Total ozone, assumed 
to be known to within ±1%, was also included as an 
additional measurement. 


Figure 2. Ozone profile retrieval using simulated 
GOME data from 255-300 nm. 

In both cases the retrieval method used was op- 
timal estimation (Rodgers, 1976, 1990) and the for- 
ward model used in the retrieval process was once 
again LOWTRAN7 (Kneizys, op. cit.). The a pri- 
ori information was assumed to be known to within 
±30% (in volume mixing ratio) of the “true" profile 
used to simulate the measurements. 

Both simulations were carried out under the fol- 
lowing conditions: nadir viewing; no clouds; no 
aerosol; over-head sun; a surface albedo of 0.3; tem- 
perature profile known exactly and spectroscopic 
data known exactly. The background atmosphere 
used was the US Standard Atmosphere (NO A A, 
1976). The results obtained under these conditions 
are discussed in the following section. 

RESULTS AND DISCUSSION 

As shown in fig. 2, the retrieval carried out us- 
ing only measurements in the 255-300 nm region, 
produces good profile information over the altitude 
range 20-60 km. In this case, height information is 
obtained by exploiting the wavelength dependence 
of effective scattering altitude for radiation at wave- 
lengths shorter than 300 nm. As noted, this tech- 
nique has been successfully applied to the SBUV in- 
strument data sets (Heath et a/., op. cit). Since the 
Huggins bands have been excluded, and all radia- 
tion at wavelengths shorter than 300 nm is scattered 
above the tropopause (Brasseur tz Solomon, 1986, 


2 



Figure 3 . Ozone profile retrieval using simulated 
GOME data from 255-340 nm, including the tem- 
perature dependent Huggins bands. 


fig. 4.3) no tropospheric profile information has been 
obtained. Conversely, when the Huggins bands are 
included in the retrieval the tropospheric profile is 
well determined (fig. 3). The relative contribution 
of tropospheric ozone to the total absorption in the 
Huggins bands may be seen in fig. 4. This shows 
the ratio of the measured backscattered spectrum to 
the incident solar spectrum, both with and without 
tropospheric ozone. It is possible to see that tropo- 
spheric ozone does make a notable contribution to 
the total absorption in the Huggins bands, thus pro- 
viding sufficient information to discriminate between 
the stratospheric and tropospheric profiles. Inclusion 
of the Huggins bands also significantly improves the 
information content of the retrievals at stratospheric 
altitudes, up to 40 km. However, these simulations 
represent the best possible situation only. In prac- 
tice, factors such as uncertainties in spectroscopic 
parameters, uncertainties in the temperature profile, 
poor a priori information, uncertainties in the radio- 
metric calibration of the instument etc. will all serve 
to degrade the quality of the profile information. 

CONCLUSION 

The potential of the ESA GOME instrument, to 
produce both stratosheric and tropospheric profiles 
of ozone, has been demonstrated by carrying out two 


Figure 4. Ratio of the measured back-scattered spec- 
trum to the incident solar spectrum , both with and 
without tropospheric ozone. 

sample retrievals using simulated GOME data. The 
simulations have shown that it is possible to obtain 
good ozone profile information from the ground up 
to 60 km by taking advantage of the temperature de- 
pendence of the ozone Huggins bands. It must be 
recognized that these simulations represent the best 
possible situation only. In practice, there will be 
many other souces of error which will increase the 
uncertainties on the retrieved profiles. In addition, 
it has been shown that using the GOME data, it is 
possible to retrieve stratospheric ozone profiles us- 
ing the traditional back-scatter ultraviolet technique 
exploited by NASA’s SBUV instruments. 
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Abstract. Cloud height and cloud coverage detection are important for total ozone 
retrieval using ultraviolet and visible scattered light. Use of the 0 2 A and B bands, 
around 761 and 687 nm, by a satellite-borne instrument of moderately high spectral 
resolution viewing in the nadir makes it possible to detect cloud top height and related 
parameters, including fractional coverage. The measured values of a satellite-borne 
spectrometer are convolutions of the instrument slit function and the atmospheric 
transmittance between cloud top and satellite. Studies here determine the optical depth 
between a satellite orbit and the Earth or cloud top height to high accuracy using 
FASCODE 3. Cloud top height and a cloud coverage parameter are determined by 
least squares fitting to calculated radiance ratios in the oxygen bands. A grid search 
method is used to search the parameter space of cloud top height and the coverage 
parameter to minimize an appropriate sum of squares of deviations. For this search, 
nonlinearity of the atmospheric transmittance (i.e., leverage based on varying amounts 
of saturation in the absorption spectrum) is important for distinguishing between cloud 
top height and fractional coverage. Using the above-mentioned method, an operational 
cloud detection algorithm which uses minimal computation time can be implemented. 


1. Introduction 

More than 50% of the Earth’s surface is covered with 
clouds. Cloud height and cloud coverage detection are 
important for a number of atmospheric measurement pur- 
poses, including total ozone retrieval and distinction be- 
tween the stratospheric and tropospheric 0 3 burdens using 
ultraviolet and visible scattered light. The Nimbus 7 total 
ozone mapping spectrometer and solar backscattered ultra- 
violet instruments use scattered solar ultraviolet light for 
ozone profile and total ozone retrieval [Fleig et al ., 1990]. 
They estimate cloud coverage from measured reflectivity at 
340 nm, assuming that the cloud top exists uniformly at a 
pressure of 0.4 atm. For improvement in environmental 
monitoring, more detailed cloud detection is required. 

To collect enough photons for the required signal-to-noise 
ratios in measurements, a spacebome spectrometer usually 
has a larger instantaneous field of view (IFOV) than radiom- 
eters like the thematic mapper (LandSat-D) and Systfcme 
Probatoire d’Observation de la Terre. The large IFOV of the 
spectrometer will often include both cloud-containing and 
cloud-free areas. High spatial resolution cloud distribution 
data are obtained by the spacebome radiometers, but it is 
difficult for them to detect cloud height. 

Figure 1 shows the calculated radiance in the 0 2 visible 
region for complete coverage in a few cloud type cases using 
the radiative transfer model MODTRAN [Berk et al ., 1989] 
in the case of 6 = 60° and <p = 0°, where 6 is the solar zenith 
angle and <p is the nadir- vie wing angle. The absolute level of 
the measured radiance depends on the cloud type. Also, it is 
difficult to distinguish cloud from Earth albedo variation. 
Thus a method providing direct information on the atmo- 
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spheric structure, such as the measurement of relative 
radiances using absorption lines, is necessary for improved 
cloud parameter detection. Use of the 0 2 A and B bands, 
around 761 and 687 nm, is such a method, one which we will 
demonstrate makes it possible to detect cloud height and 
coverage parameters simultaneously from a nadir-viewing 
instrument. Cloud altitude estimation using the 0 2 A band 
was suggested by Yamamoto and Wark [1961], and a de- 
tailed study was done by Wark and Mercer [1965]. Saiedy et 
al. [1967] determined cloud top height on a Gemini mission 
using a spacecraft-borne spectrograph camera with 0.5-nm 
resolution. Fisher and Grassl [1991] and Fisher et al. [1991] 
measured the 0 2 A band with an airborne 0.6-nm spectral 
resolution spectrometer having high spatial resolution and 
detected cloud top height with 50-m accuracy. However, it is 
rare that an entire IFOV for a satellite-based spectrometer 
(typically 100 x 100 km) is covered evenly by a cloud layer. 
Therefore cloud coverage detection is also important. We 
show that using a moderately high (approximately 4.0-cm _1 
resolution) spectrometer and measuring in several absorp- 
tion channels, both cloud top height and cloud coverage 
detection become possible. 

The study done here is a prototype for the Global Ozone 
Monitoring Experiment (GOME) and the scanning imaging 
absorption spectrometer for atmospheric chartography 
(SCIAMACHY), which are European satellite instruments 
designed primarily for atmospheric trace gas measurements. 
GOME and SCIAMACHY are both diode array-based spec- 
trometers, which measure the atmospheric spectrum with 
continuous coverage from 240 to 790 nm at a resolution of 
about 0.1 nm in the ultraviolet and 0.2 nm in the visible. 
Measurements in this region are expected to provide deter- 
minations of the atmospheric trace species 0 3 , N0 2 , and 
BrO, and the more abundant species 0 2 , (0 2 ) 2 , and H 2 0, 
and have prospects for measuring H 2 CO, S0 2 , CIO, and 


14,481 



14,482 


KUZE AND CHANCE: CLOUD PARAMETERS FROM VISIBLE O z BANDS 


| i i I | i n* |"i i i | i i i | i i i | 



600 620 640 660 680 700 720 740 760 780 800 

Wavelength (nm) 


Figure 1. Calculated radiance in the visible region, including the 0 2 A and B bands, for complete 
coverage of selected cloud types. From the top, altostratus, stratus, and cloud free. 


OCIO under special conditions. The wavelength coverage of 
both instruments includes the ultraviolet Hartley and Hug- 
gins bands and the visible Chappuis bands of 0 3 . The 
continuous spectral coverage of SCIAMACHY continues 
into the infrared, to 1700 nm, with additional detector arrays 
at 2.0 and 2.3 jim. GOME is a nadir- viewing instrument, 
while SCIAMACHY includes a nadir-viewing mode. Both 
have the capability to make measurements at the resolution 
and accuracy needed for application of the principles inves- 
tigated in the present study. GOME is scheduled for launch 
on the European Space Agency’s ERS 2 satellite in late 
1994/early 1995, while SCIAMACHY will be a part of the 
Envisat 1 payload in the late 1990s. The simultaneous 
detection of cloud parameters will help considerably in trace 
gas retrievals for these instruments. For ozone in particular, 
the concentration peaks in the stratosphere (nominally at 
—23 km), well above normal cloud altitudes. However, 
ozone in the troposphere makes a nonnegligible contribution 
to the total column density, one which may be growing as a 
result of anthropogenic activities. A major impetus for this 
study is the correction of ozone measurements in order to 
permit both better statistical evaluation of global ozone 
column measurements and determination of the climatology 
of tropospheric ozone. 

In the ultraviolet and visible regions, incident light at 
different wavelengths penetrates to almost the same distance 
inside of clouds; the effective cloud top height does not 
depend strongly on wavelength. Thus measurements of the 
0 2 A and B bands are useful for 300- to 400-nm cloud top 
height detection, appropriate to measurements in Hartley 
and Huggins bands of 0 3 , as well as in the visible Chappuis 
bands. Both GOME and SCIAMACHY will be capable of 
detecting both cloud top height and cloud coverage with 
moderate accuracy when combined with appropriate atmo- 
spheric modeling capability. The remainder of this paper will 
discuss the cloud detection capability specifically in terms of 


GOME, but the discussion is equally applicable to SCIA- 
MACHY nadir-viewing measurements. 

In this study we calculate the dependence of the measured 
spectrum on the cloud top height and a cloud coverage 
parameter (specifically, the product of the fractional cover- 
age times the cloud reflectivity, ar , divided by the scene- 
averaged albedo, y) for several cases of instrument resolu- 
tion, wavelength coverage, and cloud conditions. A 
minimization technique is applied in order to determine the 
range of values which correspond most closely to the mea- 
surement. The scene-averaged albedo can be determined 
independently by measurements out of the 0 2 bands in order 
to determine ar. Complete determination of the fractional 
coverage currently depends on independent knowledge of 
either the cloud reflectivity or the ground albedo. In future 
work we hope to be able to demonstrate techniques for 
simultaneous retrieval of cloud reflectivity and ground al- 
bedo as well as cloud top height and cloud coverage with 
instruments of sufficient resolution and sensitivity. Such 
retrieval would remove the ambiguity in the coverage pa- 
rameter and its dependence on climatological values of cloud 
top reflectivity and Earth albedo. 

2. Radiative Transfer in the 0 2 Absorption 
Region 

Figure 2 shows the 0 2 A and B spectral bands, calculated 
using the 1992 HITRAN database [Rothman et a/., 1992], 0 2 
A and B have strong absorption around 13,143 cm -1 (761 
nm) and 14,546 cm‘ l (687 nm). (Note that wavelength, in 
nanometers, is the unit of choice for measurements by the 
satellite instruments, which are dispersive spectrometers. 
Line parameter listings are normally in cm" 1 , since line 
shapes are appropriately described in a scale linear in 
energy.) A spectrometer with a spectral resolution of about 
4.0 cm -1 (e.g., GOME) will have several spectral measure- 
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Figure 2. A stick spectrum showing the 0 2 A and B band positions and intensities from the 1992 
HITRAN database. 


ments in the 0 2 bands. The radiance measured in the 
satellite geometry by spectral channel j, neglecting scatter- 
ing, is described by 


I(j) » 



fj(v)F(v) exp hi)) — 

Av 



dv 

fj(p)F( v) exp (~ST(v t 0)) — , 

Av 


(1) 


where N is the number of different cloud top height types, 
is the cloud top reflectivity for cloud type i in spectral 
channel j, r { is the coverage of type i cloud, fj(v) is the slit 
function of channel J, F(v) is the solar spectrum, i{v, h) is 
the optical depth between cloud top and satellite (propor- 
tional to total 0 2 column amount above the cloud top), h t is 
the top height of type i cloud, Ap is the full spectral width at 
half maximum of the instrument, and fy is the Earth’s 
surface diffusive albedo. Here and /3 \j are assumed to be 
constant over the small spectral size of spectrometer channel 
considered here. The optical path factor s is sec(0 4 - sec(^>). 
The cloud type definitions that we use correspond to those of 
LOWTRAN and MODTRAN [Kneizys et al 1988; Berk et 
al. y 1989). The major effect of different cloud types (except 
for cirrus, as noted below) is to change the radiance level of 
the reflected radiation. From studies such as the 
MODTRAN calculations shown in Figure I, we can assume 
that clouds are thick enough that solar light reflects at the 
cloud tops; this is the operational definition of “cloud top” 
in these studies. There are several caveats in this assump- 
tion. If a cloud layer is not sufficiently thick that determina- 
tion of a cloud top height to the — 1-km accuracy aimed for in 
this study is possible, there will obviously be an ambiguity in 
the determination. Different cloud types may have bidirec- 


tional reflectance distribution functions (BRDFs) that differ 
sufficiently that the reflectance calculations at different view- 
ing and solar zenith angles are significantly affected; this 
aspect of the problem is currently under investigation at the 
University of Bremen, Bremen, Germany, (J. Burrows, R. 
Spurr, and T. Kurosu, private communication, 1993). The 
effect of cirrus clouds on the determination of cloud param- 
eters deserves, and will get, more attention than it has 
received in this study. Our current feeling is that optically 
thin cirrus will not present a significant problem; clouds that 
are partially transmitting have the same effect when aver- 
aged over the instrumental field of view as a partial cover- 
age. Multiple Mie scattering inside of clouds may signifi- 
cantly affect the optical path; the size of this effect is 
currently being investigated at the University of Heidelberg, 
Heidelberg, Germany, with Monte Carlo multiple scattering 
calculations (H. Frank and U. Platt, private communication, 
1993). This specifically includes an investigation into poten- 
tial multiple scattering problems in cirrus. The effect of 
multiple scattering inside clouds may be the limiting factor in 
the general applicability of optical methods in cloud deter- 
minations. Radiance in the absence of absorption by 0 2 or 
other gases (e.g., out of band), without scattering, is de- 
scribed as 


IoU) = y J J F(v)fj{v) — , 

where 

-S') 


N 


r/ = 2 + p A 




( 2 ) 


For GOME, F(v) is determined from a priori information 
and on-board solar flux monitoring using a diffuser. Thus the 
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scene-averaged albedo yj out of band can be directly deter- 
mined. 

Scattering effects are considered below. In the above 
wavelength region, absorptions by gases other than 0 2 are 
small compared with the dominant 0 2 absorption; they can 
be safely neglected for cloud detection at the present level of 
accuracy. In the 0 2 A and B bands, yj is interpolated using 
the spectrum near the 0 2 A and B bands where there is no 
strong absorption. The slit function fj(v) is assumed to be 
triangular, because the spacebome spectrometers we are 
considering have the same output slit (detector) size as the 
entrance slit. In the real measurements the slit function 
must, of course, be measured with high spectral resolution. 
The vertical profiles of 0 2 number density and temperature 
are assumed to be well known. Optical depths T<y, h) 
between the satellite orbit and each height are calculated 
with high accuracy using FASCODE 3 (G. P. Anderson and 
J. H. Chetwynd, private communication, 1992). 

The measured values are convolutions of the slit function, 
/ 0 (j/), and the atmospheric transmittance. By analogy with 
Fleig et al [1990] it is convenient to describe the convolution 
in terms of a parameter Q , defined in our case as 

Q(j, hj) = J fj(v) exp (-JTfi', hi)) — . (3) 

The parameter R(v 0 , h t ), which is the ratio of satellite- 
measured radiation to solar input radiation, is then defined 
by the following, assuming F(v) is constant within a spec- 
trometer channel so that it can be taken outside of the 
integral in equations (1) and (2): 

N /O') 

R(j> h' n a if i) 


N N 

= 2 a <fiQU, h^ + j3,(i - 2 n)QU, 0) 

f- 1 1=1 

N N 

= 2 A i) + lyj - 2 “ y T i)CO'» 0) (4) 

1*1 *°i 


3. Determination of Cloud Top Height 

and Related Parameters 

3.1. Grid Search Least Squares Fit 

Cloud parameters are determined by least squares fitting 
to calculated radiance ratios. We assume there is no a priori 
information on A,* and a 1 * / r; in equation (4) and use a grid 
search method to search the parameter space of h iy a^r t to 
minimize the variance. The quantity x 1 iS defined as 

*obsO, hh ~ flcaicO> h <> 

FobsiJi ® 

where the are R values from measurement and dividing 
by the measured solar calibration, the R^c are R values for 
grid searching to determine A, and a ; r, , and M is the number 
of spectrometer channels in the 0 2 A and B band regions ( a 
and yj are assumed to be constant within each 0 2 band). 

Here x 2 not divided by the error value a 2 , because a 2 is 



M y (5) 


almost the same for each spectrometer channel. The a 2 
attained in the measurements determines the contour de- 
scribing the level of accuracy achieved in cloud parameters, 
as shown below. The ranges of r, and A* are between 0 and 
1 and between 0 and 10 km, respectively. In the present 
work, grid searching is used rather than a continuous non- 
linear least squares technique, since we are evaluating x 2 
over a space of variables with discrete values. No particular 
attempt has yet been made to optimize the technique; rather, 
we perform complete evaluations of x 2 over the range of 
coverage variable at — 0.01 resolution and over cloud-top 
height at M).l-km resolution. Because the ranges and reso- 
lutions are limited, the computation time necessary to search 
the resulting tables for minimum values is not a serious 
concern. 

For both cloud top height and cloud coverage detection, 
optical depth is an important concern. If the optical thick- 
ness is too large, it becomes impossible to detect low cloud 
accurately, because the Q values become too small. If the 
optical depth is too small, the transmittance is an almost 
linear function of a^r,- and optical depth, and it becomes 
impossible to distinguish cloud height from cloud coverage. 
If the difference of integrated optical depth between chan- 
nels is too small, it is also impossible to distinguish cloud 
height from cloud coverage. Thus for cloud parameter 
detection, both the optical depth values and contrast be- 
tween channels are important. 

Figures 3a and 3b show calculated optical depths for a 
range of typical cloud top height scenarios using FASCODE 
3: Figure 3a for the 0 2 A band and Figure 3b for the 0 2 B 
band. Around this spectral region the GOME resolution is 
about 4.0 cm” 1 . The optical depths change substantially in 
this region, so calculations must be done with high spectral 
resolution. When the spectral resolution is much wider than 
the oxygen spectral lines in the bands, the measured optical 
depth is a highly averaged or smeared-out quantity and its 
variation with detector channel is small. Therefore higher 
spectral resolution improves cloud detection. Figure 4 shows 
Q values for several A and B band channels in the case of 
Av = 4.0 cm -1 and s - 2.0. The center positions of the 
channels are 13,132, 13,144, 14,539, and 14,551 cm* 1 . 

We have calculated the variation of x 2 in the space of 
cloud top height versus the coverage parameter a x r } iy for 
representative cases of instrumental resolution, number of 
measuring channels, and cloud top heights. The results are 
given in Figure 5 as contour plots showing the logarithm of 
X 2 ; these plots are intended to give a visual sense of how well 
a minimization scheme (such as the grid search employed 
here) will succeed in the independent determination of cloud 
top height and the coverage parameter. Figure 5a shows the 
variation of x 2 for a typical GOME case, A v = 4.0 cm 1 , 
S = 1, and Af = 12, with a cloud layer at 5 km and a 
coverage parameter of atf/Vy = 0.5 (50% coverage when 
the cloud top reflectivity is equal to the Earth albedo). The 
twelve spectrometer channels chosen consist of seven 0 2 A 
and five 0 2 B band measurements. The center positions of 
the channels are 13,132, 13,136, 13,140, . . . 13,156 cm -1 for 
the 0 2 A band and 14,539, 14,543, 14,547, . . . 14,555 cm 
for the 0 2 B band. N = 1 is assumed because there will most 
often be one cloud type within the IFOV. Figure 5b com- 
pares this case with that of higher spectral resolution, Av = 
0.2 cm -1 . The center positions of the channels are 13,150.6, 
13,150.8, 13,151.0, ... 13,151.8 cm' 1 for the 0 2 A band and 
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Figure 3. Calculated nadir optical depths for the visible bands of 0 2 between space and the altitudes 0 
(top), 2, 4, 6, and 8 (bottom) km: (a) 0 2 A band and (b) 0 2 B band. 


14,555.1, 14,555.3, 14,555.5, . . . 14,555.9 cm -1 for the 0 2 B 
band. Figures 5c and 5d show cases of a lower cloud top 
height (2 km) with Ay = 4.0 cm" 1 and Ay = 0.20 cm" 1 as in 
Figures 5a and 5b. Figure 5c has several minima, but the 
proper one is significantly deeper; if the measurement accu- 
racy is very good, the solution is unique. For the lower- 
altitude cloud and lower spectral resolution there are ex- 
tended trough regions in the variation of x 2 with cloud height 
versus cloud coverage. This indicates that high accuracy is 
required in order to get a unique solution for cloud coverage. 
Generally, higher resolution is better for cloud detection, but 
there may be a tradeoff here between resolution and signal- 
to-noise ratio. However, as Figure 5c shows, cloud height is 
easily determined. For total ozone retrieval, the coverage of 


low cloud may not be critically important. From the hyper- 
bolic shapes of the contours in Figure 5c we may conclude 
that for low cloud top height and low cloud coverage, the 
cloud top height is not accurately determined. On the other 
hand, for the higher cloud top height and higher coverage 
case, both height and coverage detection accuracy is high. 
There is a basic tradeoff in cloud coverage versus cloud 
height information that stems from the fact that the most 
basic accurate information obtained by measurements in the 
0 2 bands, which is that of absorption in the optically thin 
portions of the bands, gives the average transmission or, 
equivalently, path length through the atmosphere. Any over- 
estimate of cloud coverage thus implies an underestimate in 
cloud top height. Figures 5e and 5f are moderate- and 
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Figure 4. Q values (convolutions of the instrument slit function and the atmospheric transmittance) for 
0 2 AandZ?band channels with center positions of 14,539, 14,551, 13,132, and 13,144 cm” 1 (top to bottom) 
for instrument resolution Ar = 4.0 cm -1 . 


high-resolution cases of M = 2, r x = 0.5, and h { = 5 km. 
These are calculated in order to compare results obtained 
with differing numbers of measurement channels. The reso- 
lution is 4.0 cm" 1 , and the center wavelengths of the 
channels are 13,140 and 14,547 cm -1 for Figure 5e; the 
resolution is 0.2 cm” 1 , and the center wavelengths of the 
channels are 13,151.0 and 14,555.5 cm” 1 for Figure 5f. 
Comparison of Figures 5e and 5a shows clearly that the use 
of a large number of spectral channels greatly improves the 
accuracy for moderate (i.e., GOME) resolution. However, 
for higher spectral resolution, a small number of spectral 
channels is adequate, as shown by Figure 5f. 

The above method determines ajr, . By using interpolated 
values of y Jy 0(1 - r,) may also be determined. If appro- 
priate a priori data for 0 are available, a x and r j may be 
determined individually. From the a x value it may be 
possible to estimate the cloud depth. This study is for one 
cloud type and cloud top height, but the determination of 
two or more cloud layers within the IFOV may be possible 
with the use of a grid search in more than four dimensions. 

3.2. Accuracy Requirements and Error Estimation 

The desired uniqueness of solution determines the accu- 
racy requirement. As mentioned above, for lower-height 
cloud detection, high accuracy is required. This statement 
can be quantified in a typical example. Figure 6 indicates the 
|*| variation for N = 1, M = 12, r x = 0.5, h x = 5 km, and 
kv = 4.0 cm ”* . As before, the twelve spectrometer channels 
include seven in the A band and five in the B band. The 
central contour, shown as a bold line, corresponds to the 
limit of ± 1% measurement error about the central minimum. 
This contour corresponds to ±20% error in cloud coverage 
and ±2 km in cloud top height. If the cloud coverage 
detection accuracy is required to be ±20%, the sum of 
measurement error and model calculation error must be less 


than 1%. The anticipated random errors in measurements 
are small enough to allow precise solution (the modeled 
signal- to-noise ratio for GOME is —1000), but the forward 
model calculations in the present study have uncertainties 
which may cause substantial systematic error. It is possible 
to detect the area surrounded by the contour line in Figure 6 
when the systematic errors are 1%. This area indicates the 
correlation between cloud top and cloud coverage and 
should be useful for removing cloud effects in total ozone 
determinations. To separate cloud top height and coverage 
accurately, higher spectral resolution is required. The signif- 
icant systematic error sources are discussed below. These 
items are expected to be corrected by adequate calibration. 

Measurement errors. Because R is a ratio, that of the 
Earth measurement to the solar measurement, the system- 
atic radiance errors of the instrument response largely can- 
cel. The capability for relative radiance measurements 
should be better than 1%. 

Model calculation errors include the items listed below. 
The first two items affect both y and Q value calculations. 
The third affects y, and the others affect the Q values. 

Single and multiple scattering effects. As discussed be- 
low, scattering effects must be included for accurate cloud 
determination. 

Effects of other gases. There are strong H 2 0 bands be- 
tween the 0 2 A and B bands, but they do not affect the 0 2 
A and B bands significantly. In FASCODE 3 studies we find 
the calculated optical depth due to other gases within the 0 2 
A and B band limits to be negligible. 

Scene-averaged albedo interpolation in the 0 2 A and B 
bands. The solar spectrum can be assumed to decrease 
linearly around the 0 2 A and B bands. The error caused by 
interpolation is negligibly small. 

Line parameter database. As discussed below, the line 
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Figure 5. The x 2 variation in log scale, where x 2 is the sum of squares of radiance ratio deviations, which 
is minimized in the fitting process to determine cloud-top height and coverage (see text for details). 
Estimated uncertainties are not included in x 2 \ rather, the uncertainties determine which contour is 
applicable. Results are presented for several case studies illustrating differences due to the number of 
measurement channels, AT ; the height of the cloud layer, h ; and the instrument resolution, Av. The cloud 
coverage parameter, a x r x !y = 0.5, includes a ratio of albedos, as explained in the text, (a) M = 12, 
ct\r x fy = 0.5, h = 5 km, and Av = 4.0 cm -1 ; (b) M = 12, a x r x ly - 0.5, h = 5 km, and Av = 0.2 cm -1 ; 
(c) M — 12, ajTj/y = 0.5, h - 2 km, and Av = 4.0 cm” 1 ; (d) M = 12. a x r x !y = 0.5, h - 2 km, and 
Av - 0.2 cm” 1 ; (e) M = 2, a x r 1 ly = 0.5, h = 5 km, and Av = 4.0 cm” ; and (f) M = 2, a } r x fy - 0.5, 
h - 5 km, and Av - 0.2 cm” 1 . 


intensities and pressure-broadening coefficients must be 
determined to high accuracy compared to that of the desired 
detection accuracy, or an appropriate in-flight characteriza- 
tion and calibration scheme must be used. 

Temperature profiles. According to the Work and Mercer 
[1965] absorptance calculations, the effect of temperature 
profile ambiguity on optical depth calculations in these bands 


(at low spectral resolution) is less than 1%. A calculational 
check of the temperature dependences for selected individ- 
ual line strengths confirms this conclusion for the stronger 
lines in the A band. 

Slit function. The absolute value of spectral responsivity 
is not a concern in this method. The shape of the modeled 
instrument slit function is, however, extremely important. 
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Figure 5. (continued) 


Extra absorption inside clouds. The estimated cloud top 
is always below the actual cloud top because of the extra 
absorption inside the clouds [cf. Saiedy et al . , 1967]. For our 
present purposes of correcting for ozone and other trace 
species, the important consideration is whether multiple Mie 
scattering within the cloud, of light that is eventually scat- 
tered back out into the instrument field of view, is a severe 
limitation. If the major effect of penetration into the cloud is 
replacement of the cloud top by a diffuse scattering surface 
at a lower altitude, then the technique is perfectly suited for 
such correction. This problem is not thought to be serious at 
the present (1 km) level of accuracy, but the previously 
mentioned calculations underway at the Universities of 
Bremen and Heidelberg will provide more definitive answers 
to this concern. 


3.3. Algorithm and Operation 

Before launch, the optical depths necessary to determine Q 
values are calculated for several temperature profile scenarios 
at a resolution of 0.001 cm*" 1 , and the slit function is measured 
at similarly high resolution. Look-up tables in terms of optical 
depth are then prepared at intervals of —0.1 km. 

After launch, the operational software searches for the 
minimum x 2 and estimates the cloud top heights and cover- 
age. The increments for heights and coverages are —0.1 km 
and 1%, respectively. The default number of different cloud 
types within the IFOV, Af, is 1. Thermal differences between 
ground calibration and orbit, including thermal cycling of the 
instrument through the satellite orbit, may cause the center 
wavelengths for the measurement channels to change. Thus 
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Q values must be recalculated after the on-board wavelength 
calibration. 

4. Discussion 

4.1. Cloud Reflectivity and Earth Surface Albedo 

Both a| 7-j and 0(1 - r,) are determined in this method. 
As mentioned above, in order to estimate cloud coverage 
and cloud reflectivity we will, as our initial working algo- 
rithm, employ an a priori Earth surface albedo table which 
includes both geographical and seasonal information. The 
source of the albedo database used to create this table is the 
ETOPOS data set (N OAA/National Geophysical Data Cen- 
ter, Boulder, Colorado). More sophisticated future versions 
of the algorithm should exploit the radiance data from the 


satellite instrument itself further. GOME and SC IAMACHY 
each measure the scene-averaged albedo, y J7 from 240 nm to 
790 nm, including measurements in wide spectral regions 
where atmospheric gas absorptions are small. As 0 has 
substantially stronger wavelength dependence than a } , this 
scene-averaged albedo information should help to further 
separate surface albedo and cloud coverage. Additionally, 
cloud scattering studies currently underway at the Univer- 
sity of Bremen should provide invaluable information on 
cloud top reflection coefficients coupled to cloud type and 
altitude (J. Burrows, R. Spun*, and T. Kurosu, private com- 
munication, 1993). Thus we anticipate that the proper combi- 
nation of the satellite radiance measurements and cloud scat- 
tering studies will provide the necessary leverage to separate 
albedo, cloud reflectivity, and coverage to high accuracy. 
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Figure 6. The variation in log scale; M = 12, a x r x /y = 0.5, h = 5 km, and Ap = 4.0 cm *. 


4.2. Scattering Effects 

Even though Rayleigh scattering is small in the visible 0 2 
absorption region compared with the ultraviolet region cov- 
ered by GOME, its effect in the atmospheric path must be 
considered for accurate cloud detection. Single-scattering 
radiance calculations do not depend on either the cloud top 
reflectivity or Earth surface albedo, but when multiple 
scattering is included, there is a strong dependence on both. 
Therefore to estimate scattering effects, a priori data of a 
and 0 are necessary. The input radiance 7,0'), including 
scattering effects is 

IjU) = r\ j fj{v)(a x F(v)T(s 9 v , h x ) + S(p 9 h x , a^) — 

+ (1 - r x ) ( fj(v)(pF(v)T(s, p 9 0) + SO, 0, 0)) 

J V (6) 

where T(s , v 9 h x ) is the transmittance and S(p 9 h x , x) is the 
single and multiple scattered radiance. 

In general, reflectivity of cloud tops is high in the visible 
region. Therefore scattering effects are relatively small and 
radiance calculation errors are small. Out of the 0 2 absorp- 
tion region the scattering affects yj estimation. The differ- 
ence between the radiance not considering scattering, Iq(J) 9 
and the radiance considering scattering, is less than 

1% when 0 = 0.3, because the back-scattering ratio to the 
zenith is almost equal to the Earth surface diffusive reflec- 
tivity (note that geometry must be carefully considered when 
evaluating the importance of Rayleigh scattering in non- 
GOME geometries, though). A 10% ambiguity of a priori 0 in 
the S( p 9 0, 0( p)) calculation causes an error of a few percent 
in /,(/). Therefore the y y error caused by scattering is less 
than a few percent. 

In the 0 2 absorption region, scattering affects Q value 
calculations. Because the optical path, including the 0 2 
column density, is multiplied by scattering, the effective 


optical depth calculation depends on a and 0. For low 0, 
scattering effects on Q value estimation are relatively high. 
Therefore ambiguity in 0 might cause a serious error, and a 
priori data on 0 are important for estimation of scattering 
effects. 

4.3. Spectral Band Database 

The line parameters for these bands necessary to perform 
the quantitative analysis proposed here are included in the 
HITRAN and GEISA databases, derived from measure- 
ments by Miller et aL [1969] and Giver et al . [1974]. The line 
positions are well known for both bands. The intensities for 
the 0 2 B band are well known (to better than 2%) f but for the 
0 2 A band they are known to only 4%. For both bands, 
considerable work is necessary on line-broadening parame- 
ters. Pressure broadening is the chief determinant of line 
widths and thus of radiative transfer through the bands 
below about 8 km in the Earth’s atmosphere. 

0 2 broadening of the 0 2 A band lines at room temperature 
is known to about 15%. No N 2 or air broadening studies 
have been made, and no temperature dependencies have 
been measured. The values used in the catalogs for 0 2 A 
band broadening are actually the more accurately deter- 
mined 0 2 B band width parameters (see below), a substitu- 
tion made on the presumption that the vibrational depen- 
dence of pressure-broadening coefficients is negligible. 
Studies for the magnetic dipole-allowed microwave rotation- 
less and far infrared rotational lines [Chance et aL , 1991] 
show that (1) air broadening can be substantially different 
from 0 2 broadening, (2) temperature dependence can differ 
substantially from the suggested hard sphere value of T -0 ' 5 , 
and (3) state-to-state dependence, even when sharing a 
common lower state, can vary substantially between the 
microwave rotationless and far infrared rotational lines. This 
fact, plus experience with rotational and vibrational compar- 
isons for other molecules, suggests that the 0 2 A and B band 
broadening parameters may differ at a level that is important 
for detailed geophysical studies. The cumulative error in the 
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widths under atmospheric conditions, considering all of 
these factors, may well exceed 25%. 

The 0 2 B band widths are in somewhat better shape, as 
they were performed at a later time using more fully devel- 
oped techniques, but there are still only widths (albeit good 
to 2-3%) for room temperature 0 2 broadening. Broadening 
values for this band also need to be expanded to include air 
broadening and low temperatures to be useful for realistic 
geophysical analysis. 

Q value calculations depend directly on line intensities; 
2% error in the intensities corresponds to a few percent error 
in Q value calculation. Although GOME and SCIAMACHY 
have a — 4.0-cm -1 resolution, wider than the 0 2 A and B 
band line widths, the Q values are sensitive to pressure 
broadening line widths both because the pressure broaden- 
ing affects line saturation except in the optically thin limit 
and because the slit function is triangular and the Q values 
are proportional to convolutions of the line shape and the slit 
function. The exact sensitivity depends on the line center 
location in each channel. A 3% error in line width generally 
corresponds to less than a few percent error in Q value 
calculation. 

Detailed sensitivity analyses to assess the precise effect of 
line parameter uncertainties on the cloud parameter retrieval 
accuracy are planned for the next phase of this investigation. 
This will include more detailed dependence of Q values on 
line parameters and on the measurement channel placement 
with respect to the peaks of the 0 2 lines. Simple radiative 
transfer considerations imply that the determination of the 
scene average path length is a process that depends primarily 
on the optically thin components of the observations. We 
expect that at the first level of approximation, the knowledge 
of this quantity depends linearly upon the knowledge of the 
line intensities. We also expect that in order to be able to 
overcome potentially large systematic effects due to uncer- 
tainties in line parameters, an in-flight characterization and 
calibration scheme will need to be devised in order to 
effectively utilize any data other than optically thin absorp- 
tion measurements. 

5. Conclusions 

GOME and SCIAMACHY have several channels of vary- 
ing, but moderate, optical depth in the 0 2 A and B band 
regions. We have proposed an operational method for deter- 
mination of cloud parameters that has the potential to 
produce accurate results with minimal computation time. 
This method has been primarily developed as a contribution 
to the development of the GOME and SCIAMACHY pro- 
grams. It has the potential to be effective for other satellite- 
based measurement programs as well, including geostation- 
ary meteorological satellite observations that employ larger 
IFOVs and longer integration times compared with the Sun 


synchronous satellite observations of GOME and SCIA- 
MACHY. For accurate detection, further studies of spectral 
band parameters, scattering effects, and Earth albedo stud- 
ies are required. Scattering studies currently underway 
include those of cloud-top reflection coefficients (BRDFs) 
coupled to cloud type and altitude (J. Burrows, R. Spurr, 
and T. Kurosu, private communication, 1993) and the Monte 
Carlo multiple scattering calculations of scattering inside 
clouds (H. Frank and U. Platt, private communication, 
1993). 
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Abstract. ERS2 based GOME provides UV-VIS- 
spectra of nadir scattered sunlight. The spectral obser- 
vations yielded column data of OCIO and BrO besides 
those of ozone and NO2. The observations of strato- 
spheric OCIO and BrO at solar zenith angles of 80° and 
higher in polar regions are unambiguous because the 
long optical paths give rise to strong absorptions. The 
vertical column for BrO in the equatorial region is found 
to be (3.0 ± 1.5) x 10 13 molec/cm 2 which amounts to 
16 ppt in the stratosphere. GOME data are compared 
to a number of groundbased differential optical spec- 
troscopy observations of BrO by zenith scattered sun- 
light at Sflndre Stromfjord, Greenland. Consequences 
for ozone depletion are discussed. 

Introduction 

Soon after the significance of chlorine as a sink for strato- 
spheric ozone, O3, became clear, the potential for O3 
depletion by bromine was pointed out by Watson, 1975 
and Wofsy et ah, 1975. Yung et al., 1980 and McElroy 
et al., 1986 recognized the synergistic coupling between 
bromine monoxide, BrO and chlorine monoxide, CIO, 
to destroy ozone. About half of this net reaction leads 
to ozone loss. Under conditions of high CIO concen- 
trations like during the antarctic ozone hole formation 
it may account for 15 to 30 % of the ozone loss in the 
polar stratospheric vortex (Jones et al., 1989, Anderson 
et al., 1990, Solomon et al., 1990). 

Under normal conditions the majority of inorganic chlo- 
rine in the stratosphere is tied up in relatively long 
lived reservoirs, HC1 and CIONO2, and the CIO con- 
centrations are low. Yet the photolysis lifetimes of the 
bromine reservoirs, HOBr and BrON02> are short re- 


sulting in considerable amounts of daytime BrO, a re- 
active bromine species. As a result of this partitioning, 
bromine is much more efficient than chlorine on a per 
atom basis at destroying ozone (WMO/UNEP, 1995). 
BrO in the stratosphere has been detected before by 
remote spectroscopy (Caroll et al., 1989; Roth, 1992; 
Richter et al.,1996 and references therein) as well as by 
in situ measurements (Brune et al., 1989; Toohey et al., 
1990). BrO is considered mainly to regenerate Br or 
Br2 by reactions (l)-(3). The BrO dimer (Br202) is 


less important than the CIO dimer (CI2O2)* 

BrO -I- BrO -4 Br + Br + O 2 (1) 

-4 Bri + O 2 (2) 

Br 2 + hv — ► Br + Br (3) 

Ozone is destroyed by reaction (4) and in addition the 
reaction of HO2 and BrO (5) leads to further destruc- 
tion of O3. 

Br + 0 3 -4 BrO + O2 (4) 

H0 2 + BrO -4 HOBr + 02 (5) 

HOBr 4 hv -4 Br + OH (6) 


At midlatitude the lower stratospheric ozone is con- 
sumed via this reaction at a rate comparable to or even 
exceeding that of chlorine alone (Poulet et al., 1992; 
Avallone et al., 1993; Garcia and Solomon, 1993). 

Experimental and Results 

GOME comprises four spectrometers observing the re- 
gion 240-790 nm simultaneously and having 4 diode ar- 
ray detectors. In total 3586 spectral measurements are 
made enabling the differential optical absorption tech- 
nique, DOAS, (Pemer et al., 1991) to be used for con- 
stituent retrieval. The spatial resolution is variable from 
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Figure 1. Differential absorption for BrO in the sum 
of 100 spectra at the Equator (114° W, 21° SZA) at day 
950907 , direct sun reference. The two absorption bands 
of the BrO present an example for the precision reached 
so far. The strength of the fitted laboratory BrO ab- 
sorption spectrum corresponds to 8 x 10 l3 molec/cm 2 . 

40 x 40 to 40 x 320 km 2 . Full coverage of the longitudinal 
distribution is reached in about 3 days for the largest 
swath width. GOME flies on ERS2 in a polar sun syn- 
chronous orbit with a period of about 100 minutes. 

Direct sun reference spectra are measured daily via a 
scattering plate. Figure 1 shows the differential absorp- 
tion spectrum of BrO obtained from the division of 100 
spectra at a solar zenith angle, SZA of 21° and the 
extra-terrestrial solar spectrum. Up to seven spectra 
can be fitted by a simultaneous nonlinear multiparam- 
eter procedure to the nadir spectra. In figure 1 and 
in figure 2 the residual absorption, having removed O3, 
O4, NO2, OCIO and Ring-spectrum, is plotted together 
with a BrO laboratory spectrum. The two absorption 
bands of the BrO are discernible. 

The average spectrum for about 21° SZA contains 8 x 
10 13 molec/cm 2 which amounts to a vertical column of 
(3.0 ± 1.5) x 10 13 molec/cm 2 assuming that most of the 
light is scattered in the troposphere. The differential 
absorption cross section was 7.9 x 10“ 18 cm 2 molec -1 
(223 K) at 355 nm and the spectral resolution 0.6 nm. 
A stratospheric mixing ratio of 14 ppt or 16 ppt is de- 
rived ass umin g a homogeneous distribution throughout 
the stratosphere starting from 16 or 17 km, respectively, 
assuming there is no BrO in the troposphere. 

The larger high frequency structure when using the sun 
reference in the background noise as compared with that 
for example in a spectrum which is obtained by dividing 
the sum of 100 spectra at 60° SZA at southern midlat- 
itudes by the sum of the 100 spectra at 21° SZA (Fig. 
2) is not yet fully understood but is most likely due to 
limitations in the measured Ring reference spectra. 

OCIO was also identified and measured by its bands at 
367, 376 and 385 nm (see also Burrows et al., 1996). 
For OCIO the differential absorption cross section of 


Figure 2. Differential slant column absorptions of BrO 
are considerable more precise when a sum of 100 spec- 
tra at 60* SZA at southern midlatitudes is Fraunhofer 
corrected by the sum of the 100 spectra at 21 0 SZA 
(Fig. 1) of the same orbit . The laboratory BrO ab- 
sorption spectrum shown corresponds in strength to the 
retrieved value. 

1.05 x 10- 17 cm 2 molec" 1 (223 K) at 376 nm was used. 
NO2 slant columns, SCs, were retrieved and tests have 
shown that the retrieved amounts of BrO and NO2 were 
independent of each other. Measurements around the 
globe of BrO, NO2 and OCIO (Fig.3) in September show 
high NO2 in the northern areas and lower values at 
the south pole. BrO is fairly symmetrically distributed 
while no OCIO is detected except in the vortex area 
around the south pole. Figure 4 shows the situation a 
month later. The northern NO2 shown at sunrise has al- 
ready diminished due to N2O5 formation at night while 
NO2 at the south pole has somewhat recovered. Partly 
as consequence of the now enhanced CIONO2 formation 
OCIO has decreased and in the north sporadic OCIO 
signals indicate the start of chlorine activation. BrO is 
diminished in southern hemisphere at large solar zenith 
angles most likely due to the lack of ozone in the re- 
gion where most of the bromine is situated and partly 
due to the BrON02 formation. The positive values of 
OCIO outside of the polar vortex require furter analysis. 

Validation 

From 950727 to 950808 remote sensing experiments of 
zenith scattered sun light were performed in Greenland 
by groundbased DO AS (GDOAS). BrO was evaluated 
and compared to the GOME data. Overpasses by ERS2 
at Sendre Stromfjord, SS, during that period were at 
~ 88° SZA. In Table 1 the SC densities of BrO as mea- 
sured by the nadir and by the zenith observations are 
given. The orbits of GOME with respect to SS are 
sometimes quite close and sometimes several hundred 
km to the east or west. 

On day 950728 the orbit is quite close to SS. GDOAS 
shows higher values which might be due to tropospheric 
BrO which was recognized by GDOAS unambiguously 
around noon every day until 950804. On day 950729 the 
orbit went over the ice cap to the east of SS. GOME in- 
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Figure 3. Slant column of BrO, OCIO and NO 2 and 
solar zenith angle as function of latitude of an orbit 
crossing the equator at 114 0 W at day 950907 (sun 
reference from 950905). 

dicates some variability but the agreement is fair. Day 
950730 the orbit went over the David Street close to 
Canada and GDOAS still being influenced by tropo- 
spheric BrO shows the highest values about twice as 
high as those of GOME. Next day the orbit is close to 
SS and crosses the coastal area. On this sunny and 
cloudfree day at SS both instruments show high values 
in excellent agreement indicating that tropospheric BrO 
was still present and was seen this time with GOME 
as well. On day 950801 the instruments agreed well 
again on BrO though the orbit went far away over the 
ice cap. On day 950806 the orbit went over the Davis 
Street close to Greenland. Larger slant columns agreed 
within 10-20%. In general the col umns seen by the two 
instruments agreed within error limits except that on 
occasions when the troposphere contained BrO or the 
orbit was further away from SS larger deviations might 
have be occurred. Calculations with GOMETRAN-f— b, 
a radiative transport program developed at the Uni- 
versity of Bremen, show that the differences in SC of 
GOME and GDOAS for SZA from 86 ° to 92° are less 
than 5% for a stratospheric BrO layer. An unintentional 
validation of BrO in those orbits by groundbased DOAS 
at Arrival Heights (Kreher et al.,1996) shows excellent 
agreement. In the beginning of 1996 calibration by the 
Greenland based instrument were resumed and vortex 
events above Sondre Stromfjord were studied. The eval- 
uation of the corresponding GOME spectra is underway. 


Figure 4 , Slant column of BrO, OCIO and NO 2 and 
solar zenith angle as function of latitude of an orbit 
crossing the equator at 114 0 W at day 951012 , (sun 
reference from 951012). 

Concluding remarks 

The validation of the GOME on-board ESR2, the Euro- 
pean Space Agency’s new environmental satellite launch- 
ed in April 1995, has resulted in the first space based 
observation of the atmospheric BrO. This observation of 
a key halogen species demonstrates the excellent perfor- 
mance of GOME and provides the opportunity to de- 
rive global maps of BrO. In the future studies of global, 
seasonal, interannual and longterm variations will be 
carried out. 
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Table 1 . Comparison of slant columns of BrO mea- 
sured by GOME and by groundbased DOAS (GDOAS) 
in the region of Sondre Strpmfjord , SS, Greenland. 
GOME evaluations were made with sun reference spec- 
tra taken at 950725 . For GDOAS afternoon spectra of 
the particular day were used which were normalized to 
an older reference with known content of BrO absorp- 
tion. Differentia] cross section at 355nm was 6.6 x 10" 1S 
cm -2 for GDOAS (1.1 nm resolution). 
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GOMETRAN: A radiative transfer model for the satellite 
project GOME, the plane-parallel version 
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Abstract The Global Ozone Monitoring Experiment (GOME) is a new, nadir-viewing 
instrument on the European Space Agency satellite ERS 2 which was launched in April 
1995. With diode-array detector technology, spectra of the upwelling radiance are 
measured simultaneously at thousands of wavelengths in the UV, visible, and near-IR 
regions. Inversion of these measurements by appropriate algorithms allows the retrieval 
of the distribution of gaseous and particulate constituents, which modulate the radiance 
by absorption, emission, or scattering processes. Such a retrieval algorithm requires 
an accurate radiative transfer model to describe the propagation of light through the 
atmosphere. The model GOMETRAN has been specially designed to fulfill the needs of 
GOME, but it applies generally to other downward looking space instrumentation in the 
UV, visible, or near-IR. In this paper the model is described in the plane-parallel version, 
and the optimization of computational parameters and comparisons with other radiative 


transfer models is presented. 


1. GOME Project 

The second European Remote-Sensing (ERS 2) satellite- 
was launched by the European Space Agency (ESA) in April 
1995. The scientific payload of ERS 2 is similar to that 
of ERS 1, which was launched in summer 1992 with a 
number of instruments designed to measure and assess the 
state of the Earth’s surface. However, to meet the increas- 
ingly important requirement to monitor atmospheric com- 
position and processes on a global scale, it was decided to 
add an atmospheric sensor, GOME (Global Ozone Moni- 
toring Experiment) to the already existing ERS 2 payload. 
GOME is a small-scale version of the scanning imaging ab- 
sorption spectrometer for atmospheric chartography (SCIA- 
MACHY), which has been selected for flight on the ESA 
environmental satellite ENVISAT 1 [Burrows et aL y 1988]. 
GOME was proposed to ESA in answer to a call for at- 
mospheric constituent monitoring instrumentation for ERS 

2. Its development represents a collaboration between ESA 
and the European atmospheric science community [Burrows 
etaU 1993]. 

The main objective of GOME is the global assessment of 
the state of the atmosphere’s composition. In recent years 
it has been recognized that both natural (e.g., volcanic erup- 
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tions) and anthropogenic activities (e.g., industrial pollution) 
may significantly influence the behavior of the atmosphere. 
The most widely known effects are the Antarctic spring- 
time ozone hole phenomenon, tropospheric pollution in and 
around industrialised areas, and the potential global climate 
warming. While ground-based measurements are usually 
sufficient to monitor well-mixed species with relatively large 
lifetimes (such as O2 and CO2), it is necessary to use satel- 
lite instrumentation to obtain a global picture of short-lived 
constituents (such as O3 and NO2), which are characterized 
by a higher horizontal and vertical variability. 

GOME is a nadir-viewing spectrometer which measures 
solar radiation backscattered from the atmosphere and re- 
flected from the Earth’s surface. The ERS 2 satellite fol- 
lows a Sun-synchronous 100-min polar orbit. The spectrum 
is measured in four separate bands from 240 to 790 nm at 
more than 3500 spectral channels, with intervals of 0.1 1 and 
0.2 nm (see Table 1). GOME is building on the expe- 
riences with the total ozone mapping system (TOMS) and 
the solar backscatter ultraviolet (SBUV) instrument on the 
Nimbus 7 satellite [Heath et al , 1975], which detected the 
solar backscattered radiation in 4 and 1 1 channels, respec- 
tively, with a width of 1 nm each. Because of the limitation 
of a small number of channels, all in the ultraviolet wave- 
length range, TOMS and SBUV constituent retrievals were 
confined to column amounts and profiles, respectively, of 
ozone only. Although the primary objective of GOME is, 
as its name suggests, also the retrieval of ozone amounts, 
the specific characteristics of the instrument enable informa- 
tion about a number of other atmospheric constituents to be 
identified in the spectrum. The gases that are potential can- 
didates for a retrieval include O3, O2, O4, the nitrogen com- 
pounds NO, N0 2 , and NO3, the halogen species OCIO, CIO, 
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Table 1. GOME Instrument 


Wavelength coverage 240 to 790 run in four bands l t 240-3 16 nm; 2, 311-405 nm; 

3, 405-61 1 nm; 4, 595-793 nm; 

Spectral intervals 0. 1 1 nm (bands 1 and 2), 0.20 nm (bands 3 and 4) 

Spectral resolution 0.22 nm (bands 1 and 2), 0.3 1 nm (bands 3 and 4) 

Instantaneous field of view 0.143° across track, 2.860° along track 

Flight altitude 780 km 

Ground speed 7 km/s 

Scan angle variable, with a maximum of ±31° 

Scanning period 4.5 s (forward scan), 1.5 s (fly back scan) 

Readout time 12 s below 312 nm, 1.5 s above 312 nm (default values) 

Spatial resolution on ground 320 x 40 km 2 (at 1.5 s readout and ±31° swath) 

960 x 80 km 2 (at 30 s readout and ±3 1° swath) 


and BrO, plus HCHO, SO2, and H2O, some of which might 
be detectable only under perturbed atmospheric conditions 
as, for example, ozone hole formation, pollution episodes, 
etc. In addition, the distribution of particulate substances 
(aerosols, tropospheric clouds, polar stratospheric clouds) 
and Earth surface reflectance properties will be investigated. 
While it will be possible to retrieve vertical profiles of O3 
and possibly of NO2 and aerosols, the retrieval of other 
gaseous constituents will most likely be limited to column 
amounts, i.e., concentrations integrated vertically over the 
atmosphere. 

The simultaneous retrieval of a multitude of parameters, 
which may be partially cross-correlated with each other, is a 
challenging task and requires a sophisticated retrieval algo- 
rithm. The retrieval method applicable to this task and cho- 
sen for GOME is based on “optimal estimation,” in which 
the measurements are supplemented by a priori information 
about the parameters to be retrieved, in order to achieve a sta- 
ble and unique inverse solution [Rodgers, 1976]. (It should 
be mentioned here that the first operational GOME retrieval 
algorithm is based on a different method and will be re- 
stricted to the retrieval of total column amounts of O3 only. 
At present, the application of the optimal estimation tech- 
nique for GOME is at a scientific level but will be considered 
for the second generation of the operational data analysis.) 
Optimal estimation is an iterative process in which the atmo- 
spheric retrieval parameters are adjusted to obtain the best 
fit between a simulated radiance spectrum and the measured 
spectrum of the radiance. The calculation of simulated ra- 
diances requires a radiative transfer model which describes 
the propagation of solar radiation through the atmosphere by 
taking into account scattering, absorption, and emission pro- 
cesses in the atmosphere, and the reflection of the Earth’s 
surface. 

The development of a GOME-specific radiative transfer 
model, GOMETRAN, is the subject of this paper. The re- 
quirements for such a model will be discussed in section 2, 
which is followed by a description of the general concept of 
GOMETRAN in section 3, and an exposition of the math^ 
ematical principles in section 4. The optimization of com- 


putational parameters which determine accuracy and speed 
is described in section 5 for the plane-parallel version of 
GOMETRAN. Section 6 presents comparisons with other 
plane-parallel radiative transfer models. Section 7 concludes 
with a summary. 

2. Requirements for a GOME-Specific 
Radiative Transfer Model 

Any radiative transfer model in an “optimal estimation” 
retrieval algorithm must calculate two kinds of output: first 
a simulated spectrum of the radiance as received at the satel- 
lite, and secondly the so-called “weighting functions.” The 
latter are the derivatives of the radiance with respect to the 
parameters to be retrieved (e.g., concentrations of O3 at vari- 
ous atmospheric levels, column amounts of gases, etc.). Ex- 
isting radiative transfer models are not set up for the con- 
venient calculation of weighting functions (this is because 
of the inherent mathematical structure). With these models 
the only way to derive weighting functions is to perturb the 
retrieval parameters singly and calculate the resulting radi- 
ance; the difference to the radiance simulated without per- 
turbation yields the desired weighting function. Depending 
on the number of parameters to be retrieved, the radiative 
transfer model has to be executed many times for each itera- 
tion step of the retrieval. 

The efficient calculation of weighting functions is one of 
the most important (but not the only) criterion for the selec- 
tion of a model applicable to GOME. The ultimate version 
of such a model must contain the following features: ( I ) 
spectral coverage of the full GOME range between 240 and 
790 nm at intervals corresponding to the instrument detec- 
tor pixels, namely 0. 1 1 nm in the ultraviolet and 0.20 nm in 
the visible; (2) inclusion of all relevant atmospheric gases 
with absorption structures in the GOME spectral range, i.e., 

0 3 , NO2, NO3, CIO, OCIO, BrO, H 2 0, O2, 0 4 , HCHO, and 
S0 2 ; (3) inclusion of all relevant emissions features, e.g., 
NO y bands; (4) inclusion of all atmospheric constituents 
with scattering properties in the GOME spectral range, i.e., 
of molecular (Rayleigh) scatterers and particulate scatter- 
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ers (aerosols and clouds); (5) an exact treatment of multiple 
scattering in the mathematical solution of the basic radiative 
transfer equation (RTE) in order to achieve accurate simula- 
tions of radiance and weighting functions; (6) an adequate 
description of the land surface spectral reflectance and of 
the bidirectional reflectance over the oceans; (7) spherical 
and refractive geometry to ensure global applicability of the 
model; and (8) high computational speed for the calculation 
of radiance and especially weighting functions; this is a re- 
quirement for future real-time data processing algorithms. A 
review of the literature showed that none of the existing and 
publicly available models could fulfill these requirements. 
The models LOWTRAN and MODTRAN both involve a 
simplified approach to describe multiple scattering and are 
therefore not accurate enough [Isaacs et al. , 1989, 1987; 
Kneizys etal> 1986; Berk et al , 1989]. FASCODE is a high- 
resolution model designed to take into account the line char- 
acter of gas absorption above 560 nm [Ridgway et al , 1982; 
Isaacs et al , 1989, 1987; Clough et al, 1986, 1989]. The 
calculation of radiative transfer by a line-by-line approach 
is too time consuming for GOME. The model DISORT has 
acceptable accuracy and speed provided that only the cal- 
culation of the spectral radiance is considered [Stamnes and 
Swanson , 1981; Stamnes and Dale, 1982; Stamnes and Con- 
klin , 1984; Stamnes et al, 1988]. 

3. General Concept of GOMETRAN 

The decision to develop a GOME- specific model, GOME- 
TRAN, was driven by the necessity to calculate efficiently 
the weighting functions in addition to the radiance spectrum. 
This requirement resulted in the selection of the finite differ- 
ence method as most suitable for this purpose [. Barkstrom , 
1976; Lenoble , 1985]. This method is amenable to a lin- 
earized description of the change in radiance caused by vari- 
ations in any of the atmospheric parameters in the radiative 
transfer equation , i.e., the weighting functions. This pro- 
vides an efficient way to perform sensitivity studies or to use 
the model in any retrieval algorithm which requires weight- 
ing functions for the fitting process (least squares, optimal 
estimation, etc.). 

GOMETRAN is available both in a plane-parallel and in a 
spherical version, the latter including refraction. The plane- 
parallel model is applicable to solar zenith angles below 75° 
and hence must exclude polar regions from its geographi- 
cal coverage. This model has been optimized and validated 
and is presented in this paper. The spherical version is still 
undergoing a thorough investigation; it is actually a quasi- 
spherical solution, whereby only the solar flux entering the 
RTE is calculated in a spherical geometry, while expressions 
for single and multiple scattering correspond to those in the 
plane-parallel RTE (see section 4). It is anticipated that this 
approach is suitable for solar zenith angles below 93° (which 
is the largest possible angle occurring during a GOME polar 
orbit). 

With regard to the wavelength coverage and spectral reso- 
lution, GOMETRAN is a generic model and can therefore be 
matched to the design of the GOME instrument. The model 


can be operated either in a monochromatic mode, where at 
each wavelength a time-consuming solution of a system of 
equations has to be performed, or in a polychromatic mode, 
where radiance calculations for a larger number of wave- 
lengths are based on a single evaluation of this solution. This 
second mode of operation is less accurate but computation- 
ally efficient (see section 5.5). 

With respect to the mathematical implementation of ab- 
sorption and emission features, GOMETRAN has been set 
up as a band model. This implies that line absorption spectra 
for gases such as O2 and H2O need to be preprocessed to a 
form suitable as input for GOMETRAN. Absorption spectra 
for other GOME-relevant gases can be taken directly from 
laboratory-measured cross sections. This refers to O3, O4, 
N0 2 , NO3, CIO, OCIO, BrO, HCHO, and S0 2 . 

Multiple scattering by molecular and particulate substan- 
ces can, in principle, be treated accurately by the finite differ- 
ence method. There are, however, two features in the current 
version of GOMETRAN that are worth mentioning. The 
first is the description of the aerosol scattering phase func- 
tions by the Henyey-Greenstein formula [ Lenoble , 1985]. 
The Henyey-Greenstein factor in this formula determines, 
in conjunction with the scattering angle, unambiguously the 
value of the phase function. Highly anisotropic phase func- 
tions cannot be characterized accurately with this approx- 
imation. The option to read phase functions from user- 
provided databases will be included shortly. The second note 
refers to clouds in GOMETRAN. At present, clouds have to 
be treated as particle layers of high optical thickness. In- 
side the cloud the radiance I changes rapidly with penetra- 
tion depth z, with correspondingly high values for the gra- 
dient dl/dz. In the finite difference method this gradient 
is approximated by a finite difference AI/Az between two 
discrete atmospheric levels used in the model (see also sec- 
tion 4). To follow radiative transport through the cloud accu- 
rately, a very large number of atmospheric levels is required 
and this makes the model too slow for any real-time data 
processing. An alternative approach is currently under in- 
vestigation for optically thick water droplet clouds, whereby 
the cloud top is treated as a bidirectional reflecting surface 
(71 Kurosu , private communication, 1994). The Earth’s sur- 
face can be selected to be either a Lambertian, wavelength- 
-independent reflector or to be described by a bidirectional 
surface reflectance. 

4. Mathematical Principles of GOMETRAN 

4.1. Calculation of Radiance 

The basis for any radiative transfer model is the integro- 
differential radiative transfer equation for the radiation field 
I. The parameters that are needed to describe the RTE are 

(1) the atmospheric composition, given by vertical concen- 
tration profiles of the gaseous and particulate constituents; 

(2) the Earth’s surface spectral reflectance; (3) the viewing 
geometry, defined by the solar zenith angle, the line-of-sight 
angle (which is closely related to the intrument scan angle), 
and the relative azimuth angle between the Sun and the line- 
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of-sight direction at top of atmosphere; and (4) a database of 
spectral parameters, including the absorption, emission, and 
scattering properties of gases and particulates. In a plane- 
-parallel atmosphere the radiation / depends only on three 
spatial coordinates: the altitude z in the atmosphere (an al- 
ternative is the optical depth), the zenith angle 0, and the 
azimuthal angle <p. The solar zenith angle Go is a parameter 
in this case. The RTE can now be expressed in the follow- 
ing form, where p is the cosine of the zenith angle 0 and 
where the dependence on the wavelength X is omitted. For 
simplicity a Lambertian surface albedo is assumed 

M^sl = _ a 

2n i (1) 

0 -1 

where 

c(z) =a(z) + b(z) 

a(z) = a M (z) + X a ' n i(z) ( 2 ) 

i=i 

b(z) = b M {z) + b R (z) 

and the boundary conditions at the top and the bottom of the 
atmosphere lead to 

I~ (zo,ft, <p) = 5(/i - p 0 , 9 - <P 

2 it 1 

/ + ( 0,/i, <p) = \ f I r (0,//, <p W<V 

0 0 

where /~(zo>/*,tp) is the downwelling radiance at zo and 
7 + (0,ji,<p) is the upwelling radiance at the ground. The other 
symbols used in the equations have the following meaning: 

/ total radiance (i.e., radiation field); 
z altitude; 

zo altitude at top of atmosphere; 
p } pi cosines of zenith angles 0, & ; 

piQ cosine of the solar zenith angle Go; 
q>, cp 7 azimuthal angles in relation to the line-of-sight 
projection on the Earth’s surface; 

(po azimuthal angle of the Sun in relation 
to the line-of-sight projection 
on the Earth’s surface; 
k number of trace gases; 
c total extinction coefficient; 
a total absorption coefficient (sum of trace gas 
and particle absorption coefficients); 
bR Rayleigh-scattering coefficient; 
bM particle-scattering coefficient; 
b total scattering coefficient (sum of Rayleigh- 
and particle-scattering coefficients); 

Pr Rayleigh-scattering phase function; 

PM particle-scattering phase function; 
p total scattering phase function (weighted 


sum of Rayleigh- and particle-scattering phase 
functions, i.e. (b R p R + b M p M )l{b R + b M ))\ 

( Xi trace gas absorption cross sections 
for all relevant gases; 
rii concentrations of trace gases; 
um particle absorption coefficient; 

itFpo solar flux at the top of the atmosphere; 

A spectral reflectance (Lambertian) of 
the Earth's surface. 

Various mathematical techniques have been developed to 
obtain a numerical solution for the radiance I from the RTE. 
Each method involves approximations which determine the 
accuracy and speed of the numerical algorithm. The “fi- 
nite difference method” chosen for GOMETRAN combines 
high accuracy with reasonable computational speed [ Leno - 
ble f 1985]. 

This method is based on an expansion of the scatter- 
ing phase function p( y) in terms of Legendre polynomials, 
where y is the scattering angle between the directions of in- 
cident and scattered radiance. (It should be noted here that 
the expansion with Legendre polynomials is not specific to 
the finite difference method but may be used in association 
with a number of methods to solve the RTE.) Using the ge- 
ometrical relation between cosine y and the cosines of the 
zenith angles p and pi and the azimuth angles (p and (p 7 for 
the directions of incidence and scattering, the phase func- 
tion can be expanded as a Fourier series of the cosine of the 
azimuthal angle 

p(z,p,p ,<?,<?') = p°(z,p,p) 

M (4) 

+2 £ /7 m (z,//,^')cos(OT(<p-(p')) 

m= 1 

where the harmonic moments p m (z y p,pl) can be described 
in terms of associated Legendre polynomials. The radiance 
I can be expanded correspondingly into a Fourier series of 
the cosine of (<p — (po). 

The combination of the two Fourier series for the phase 
function and the radiance in conjunction with the applica- 
tion of the orthogonality property of the integrals over the 
azimuthal angle <p in the RTE results in an elimination of the 
dependence on <p. The RTE reduces to a set of equations 
containing the Fourier components of the diffuse upwelling 
and downwelling radiances and Qg the subscripts m 
and diff are omitted for all the following equations for rea- 
sons of simplicity: 

= -c(z)I + (z,p) 

+ J p(z,/J,p')I + {z,p')dp' 

0 

+ -j - J p(z,p, -p)r(z,p')dp 
0 


i 

is 
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b(z) 

+ - 7 -p(z,fi,-tto)Fe «> 

4 

+ ^AFp 0 e~% J p^p^e'*^ dp' 


P ~~j~ +c(z)H(z,ft) 


/i^^ = -c(z ) r(z,^ ) 


+ j p^-p^) I+ ( z ^) d p' 

0 

+ > ^ 2 ~ / 


biz) , _l& 

+ —f-p{z,-p,-po)Fe *> 


J P~(z,P,p')H{z,p')dp' 
0 

b(z) _ — 1 

= ~£- Fe m 2 P ( Z 'P'P°) 

+ b -^AFp»e~% 

P~ {.Z,p,p')e~^ 1 dp 

0 

= Q~(z,p) 


with the abbreviations 


P + {z,p,f/) = p{z,pj) + p{z,p, -p) 
P~{z,p,p') = p{z,p,p') - p{z,p, -p') 


b(z) f 

-I — -^-AFpoe j p(z,-p,p')e * dp The boundary conditions now correspond to 


where T is the vertical optical depth and to is the optical 
depth of the entire atmosphere. For higher-order Fourier 
components (m > 0 ) the surface reflection term in this equa- 
tion is zero. The boundary conditions at the top and at the 
bottom of the atmosphere yield 


J(zo,p)-H{zo,p) = 0 


for m > 0 


l ( zo,p ) = 0 


for m > 0 


1 

J{0,p) -2 aJ J(0,p')p'dp' + H (0 ,p) for m _ 0 
0 

1 

-2 aJ H(G,p')p'dp = 0 


1 

I + (Q,p) = 2A J r(0,p')p dp for m = 0 ( 6 ) 


J(0,p) + H(0,p) = 0 


for m > 1 


l + (0,p) = 0 


for m > 1 


The next step of the finite difference method is to introduce 
two new variables, J and H : 

J(z,p) = \(I + (z,p) + I~(z,p)) 

z (7) 

H{z,p) = l+ (z,P ) - r(z,p)) 

and equation (5) now reduces to 


In the next steps discrete coordinates for the altitude z and 
the cosine of the zenith angle fi are selected and the integrals 
over the zenith angle fi in the RTE are replaced by Gaus- 
sian quadrature summation over the “discrete” streams. (In 
conjunction with the expansion of the phase function into 
a series of Legendre polynomials, the Gaussian integration 
method minimizes the residual error.) As a result, the gradi- 
ents dH/dz and dJ jdz are replaced by finite differences be- 
tween two discrete altitude levels, A H / A z and AJ/Az . These 
procedures are described in detail by [Lenoble> 1985]. Thus 
the RTE is obtained in a matrix form: 


djn^L + M 
dz 

1 

-hQ J p + (z,p,f/)J{zJW 
0 

= Fe~ l-p + (z,p,po) 

4 l 

+ ^-AFp oe~«' 

^ J P + {z,P,p)e~ X '^ dp' 

0 

= Q + (z,p) 


CL = R (11) 

All elements of the matrix C and of the vector R are known 
functions of the atmospheric, surface, viewing, and spectro- 
scopic parameters, while the vector L on the left side con- 
tains the unknown quantities H and J for each altitude level 
and zenith angle direction. The solution can be obtained by 
decomposition of the matrix C into lower and upper trian- 
gular matrices. This so called LU decomposition is required 
for each of the Fourier components. To obtain the radiance 
received at the satellite, all Fourier components have to be 
added up. The RTE considers only the “diffuse” part of the 
radiation, which has been scattered at least once. The direct 
contribution by the solar flux being reflected at the Earth s 
( 8 ) surface into the direction of the satellite is added separately. 
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The matrix C dimensions are of the order of N x N, where 
N is the product of the number of discrete atmospheric lev- 
els plus 2 with the number of discrete zenith angles (i.e., 
streams). Typical values of matrix dimension is of about 400 
x 400. To save computational time, a procedure is adopted 
to avoid matrix LU decomposition at every single wave- 
length. The need for a new decomposition when stepping 
through the wavelength range is judged by the change in ma- 
trix norms between the previous and the new wavelength. If 
the change is below a certain limit, the previously decom- 
posed matrix is used in conjunction with an appropriate cor- 
rection procedure, which is based on first-order perturbation 
theory. 

4.2. Calculation of Weighting Functions 

For the derivation of weighting functions it is necessary 
to express the relation between changes of the atmospheric 
and surface parameters x, b , c and A on the one hand and 
the radiances H and J on the other hand. For this purpose, 
new perturbed atmospheric and surface parameters are de- 
fined as p f = p + 5p, where p represents x, b , c, or A. The 
correspondingly perturbed J and H function can be written 
as 

H'(z>p,X) = H{z y p y X) + SH{z>p,X) ^ 

J'(z,p y X)=J{z y p,X) + §J(z y p,X) 

Using the perturbed values H\f ,c* ,b', T / , and A' in equation 
( 8 ) and in boundary condition ( 1 0) instead of H, J, c, b, X, and 
A and ignoring second-order terms, the first part of equation 
( 8 ) can be written as 

p 7 — - + c(z)5/(z,/i) - 

dz 

l 

[ p + {z y p y p)&{z y p)dp = 

0 (13) 

-8c(z)7(z,/i) + 

1 

J p^(z y pMz y pW + ^(z^ 

0 

the variation of the source term Q+ can be expressed in the 
following way: 

8Q + (z,p) = ^-Fe~^ l -p + (z,wo) 

-^-Fe~^ ^-p + {z,p,po)8x 
4 2po 

+^-AF m e-% J VV 

o 

-^pAF<r« 8 to J p + {z,p,p)e^^ dp/ (14) 


+ ~^-AFpoe~ « (5r - 6 to) / ^ 

0 

+ ^^-AFpoe~^8A J p + (z,p,p)e ^ dp 

o 

where the following expansion has been used: 


e 


-Si 1 

e e *>5xo 

M> 


Exchanging J and H and replacing by p~ and Q + by 
Q~ t a similar equation results for the differential dbJ/dz - 
The boundary conditions are now given by 


&J{zo,p)-5H {zo,p) = 0 

i 

bJM-iAjbJMp'df/ 

0 

+ 5H(0,p) + 2A J M{0,p')p'dp' = 
0 
l 

2&4 J (J(0,p')-H{0,p))p'dp 
o 

S/(zo,/i) + 5//(zo,f) = 0 


for m > 0 


for m = 0 


for m > 1 


The variation of Sx, 5b, and 5c arises from changes in aerosol 
scattering and absorption coefficient, Rayleigh-scattering co- 
efficient, and trace gas absorption coefficients. These param- 
eters are dependent on wavelength and are not suitable as 
parameters for weighting function calculation. New param- 
eters, which are independent of wavelength, are required: ( 1 ) 
for trace gases instead of absorption coefficients the vertical 
profiles of the trace gas concentration n[(z) have been cho- 
sen; ( 2 ) for aerosol the vertical profile of absorption and scat- 
tering coefficients at three specified reference wavelengths 
have been selected and linear interpolation has been used for 
the other wavelengths: 

for X[ <X < X 2 : 

b\f{z y X) = (1 - Yi)bA/(z,^i) + YibA/(z,^ 2 ) 
for X 2 <lX. < A 3 : (17) 

bM(z,X) = (1 - yi)bM(z,X 2 ) + Yz^Af (2,^3) 

with Yi = (A. — Xi)/(X .2 — ), Y 2 = (X- X 2 )/(Xi- X 2 ) and 

the same equation for aM(z,A)‘> (3) for Rayleigh scattering 
the wavelength-independent part of the scattering coefficient 
is selected: 

bR(z y X) = ^ R (X)r(z) (18) 

where r(z) is a function of temperature and pressure. 

Thus all atmospheric parameters are now described by 
ni(z) (i=l, k) t bM(z y X\) f bM{z y X 2 ),bM{z,X3),aM{z,Xi), 
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om(zM), a M {z,h), r(z), and A. It should be noted that the 
number of reference wavelengths h and, correspondingly, 
the number of aerosol parameters can be readily extended if 
considered necessary. 

As different atmospheric parameters have different orders of 
magnitude, a normalization is used which is given here for 
the example of trace gas concentrations: 

«Wt> = ^ i= 1, k (19) 

Using the normalized parameters, equation (13) takes the 
form 

^ dhH{z ,j i) + c(z)&/(2 ^) _ j p+izM'WizJW 

o 

P (20) 

= X v+(z',z,nWj(z')<h' 

'= ] 0 


as an example the functions for dimensionless 

variation of the vertical profile of the trace gas concentration 
is as follows: 


vt(z',z,Ar) = -a.j(z' ,tynj(z') |y(z',p)8(z - z' 

-I -p + (z,fJ,no)d{z,z) 

4 2fio 


b(z ) -Si 
4 


j p + {z,fJ,v')e ^ d P ( 21 ) 
0 

-»/ (1 } 


where P is the number of parameters; 8(z - z') is the Dirac 
delta function, and d(z,z ') is the Heaviside function: 


d(z,z') = 0 0 < z' < z 
d(z,z') = 1 z < z' < zo 


In the following steps the integrals in the above equation are 
replaced by sums with suitable Gaussian weights, and the 
differentials dSH/dz and dhJjdz are replaced by the differ- 
ences A6///Az and A8//Az. This is analogous to the treat- 
ment of the basic radiative transfer equation in section 4. The 
previous equation can be expressed now in matrix form: 

C8L = B8| (23) 


where vector 8^ consists of all vectors 8 which are rep- 
resented the functions 8^(z) in equation (20) for the finite 
altitude grid, the vector 8 L contains in an alternating way el- 
ements of 8// and 87, and the matrix C is identical to the one 
in equation (11). 

Equation (23) describes the linear relationship between 
variation of radiance and the variation of selected atmo- 


spheric parameters. This equation is, therefore, the linear 
finite difference radiative transfer equation. The solution of 
equation (23) is 

8l = C-‘B8| (24) 

As is evident from equation (24) and the above definition 
of the weighting functions as the derivatives of the radi- 
ance with respect to the atmospheric parameters, the ma- 
trix product CT l B is the finite difference approximation of 
these derivatives. It should be noted that it is not necessary 
to perform a new decomposition of matrix C to obtain the 
weighting functions but only to calculate the matrix B. A 
detailed discussion of the weighting functions generated by 
GOMETRAN is given elsewhere (V.V. Rozanov, A linear ra- 
diative transfer model for the interpretation of measurements 
of scattered and reflected solar radiation of the Earth’s atmo- 
sphere, submitted to Journal of Quantitative Spectroscopy 
and Radiative Transfer, 1996). 


5. Optimization of Computational Parameters 

The “finite difference method” chosen for GOMETRAN 
represents a numerical solution of the radiative transfer equa- 
tion (RTE). The accuracy of the solution is limited by the 
need to adopt discrete values for the height and angular co- 
ordinates and the expansion of variables in series of finite 
length. The values assigned to the related computational pa- 
rameters determine the accuracy. For GOMETRAN the ob- 
jective is an accuracy of 1 % for the simulated radiance. The 
ensuing requirements on the computational parameters have 
been investigated for the plane-parallel version of GOME- 
TRAN and are discussed in the following. 


5.1. Number of Fourier Terms 

In GOMETRAN the radiance is expanded as a Fourier se- 
ries in the azimuthal coordinate. The terms in this series are 
calculated consecutively, so that an automatic convergence 
criterion can be applied; this criterion determines the num- 
ber of Fourier terms required to achieve a given accuracy. 
An example for the convergence of radiance at five different 
wavelengths in the GOME spectral range is shown in Figure 
1 . Four Fourier terms are sufficient to achieve an accuracy 
of better than 1% in the visible wavelength range. An even 
lower number of Fourier terms is required in the ultraviolet. 
The reason for this is that the backscattered radiance in the 
ultraviolet originates from the upper parts of the atmosphere 
where molecular scattering, which has a comparably low de- 
gree of anisotropy, is dominant. This results in a low angular 
dependency in the radiance. In visible wavelength regions, 
where aerosols contribute more strongly to the radiance, the 
more anisotropic aerosol scattering phase function leads to 
higher anisotropy in the radiance and hence a larger number 
of Fourier terms is required. 


5.2. Number of Legendre Polynomials 

A second important computational parameter related to 
the azimuthal coordinate in the RTE is the number of Leg- 
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Figure 1. Dependence of the radiance on the number of 
Fourier terms. The simulated radiance is displayed at five 
wavelengths as a function of the number of Fourier terms 
used in the simulation; the radiance is normalized to that ob- 
tained using eight Fourier terms to illustrate the degree of 
convergence with respect to this “quasi-infinite” value. At- 
mospheric conditions: latitude 45°N, time of year July, with 
O 3 , N0 2 , temperature, and pressure calculated from the MPI 
two-dimensional atmospheric model [B92]; four different 
aerosol types: maritime between 0 and 2 km, tropospheric 
between 2 and 10 km, background stratospheric between 10 
and 30 km, and mesospheric between 30 and 60 km; sur- 
face albedo 0.3. Conditions for the viewing geometry: solar 
zenith angle 75°; relative azimuth angle between the Sun 
and the satellite 60° (as defined by the line of sight); line- 
-of-sight zenith angle 30°. Conditions for the computation: 
49 atmospheric layers, 7 stream zenith angles, 20 Legendre 
polynomials. 



endre polynomials. This number determines the accuracy of 
the description of the scattering phase function as series in 
Legendre polynomials. As with the number of the Fourier 
terms, there are significant differences between the behavior 
in the ultraviolet and in the visible (for an example, see Fig- 
ure 2). The molecular (Rayleigh) scattering phase function 
is fully described by the sum of a zeroth-and second-order 
Legendre polynomial. Thus the ultraviolet radiance, which 
is dominated by molecular scattering, remains almost un- 
changed for scattering expansions with three or more Legen- 
dre polynomials. In the visible the anisotropy of the aerosol 
phase function requires a higher number of Legendre expan- 
sion terms. For the aerosol types studied so far (maritime, 
tropospheric background, and stratospheric background) an 
expansion with 15 to 20 terms is sufficient for the 1% accu- 
racy objective. 

However, these results for aerosols were obtained using 
the relatively smooth Henyey-Greenstein phase functions, 
and more realistic particulate phase functions (especially 
those for water or ice particulates) would require many more 
terms in the Legendre expansion. Sharp forward peaks in 
the phase functions can be treated separately after an ap- 
proximation by a delta function, and this helps to lower the 
Legendre number. 

5.3. Number of Zenith Stream Angles 

The numerical solution of the RTE in GOMETRAN re- * 
quires a discretization of the coordinate of the zenith angle. 



Figure 2. Dependence of the radiance on the number of 
Legendre terms. The simulated radiance is displayed at five 
wavelengths as a function of the number of Legendre terms 
used in the simulation; the radiance is normalized to that ob- 
tained using 20 Legendre terms to illustrate the degree of 
convergence with respect to this “quasi-infinite” value. Con- 
ditions for the atmosphere, the viewing geometry, and the 
computation are as in Figure 1 but with four Fourier terms 
and the number of Legendre terms being varied. 


The “streams” of radiation propagating through the atmo- 
sphere are calculated in these discrete directions. The direc- 
tions are then the abscissae used in the Gaussian quadrature 
of the zenith angle integrals in the RTE, which are replaced 
by summations. The number of streams should be at least 
half the number of Legendre polynomials to minimize the 
error caused by the Gaussian quadrature, which replaces the 
integral over the zenith angle in the RTE by a summation 
over discrete streams. The number of streams has a signif- 
icant effect on the computational time because the dimen- 
sions of the matrix, which has to undergo the LU decompo- 
sition, are proportional to the number of stream angles. 

Five different wavelengths were selected to study the num- 
ber of streams required to achieve a 1% convergence. Fig- 



Figure 3. Dependence of the radiance on the number of 
zenith stream angles. The simulated radiance is displayed 
at five wavelengths as a function of the number of stream 
zenith angles terms used in the simulation; the radiance is 
normalized to that obtained using nine stream angles to illus- 
trate the degree of convergence with respect to this “quasi- 
infinite” value. Conditions for the atmosphere, the viewing 
geometry, and the computation are as in Figure 1 but with 
four Fourier terms and the number of stream angles being 
varied. 
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ure 3 shows that seven streams (for the half space) are suffi- 
cient in the visible range, whereas in the ultraviolet the num- 
ber of streams is not critical due to the dominance of single 
scattering in this spectral range. Therefore in the ultraviolet 
the multiple-scattering and the surface reflection terms in the 
RTE can be neglected, thus making the Gaussian quadrature 
of the integrals redundant. The remaining single-scattering 
part of the RTE can be solved directly for a stream into the 
line-of-sight direction and is hence not dependent on the 
zenith stream angles. 

5.4. Selection of Altitude Levels 

The selection of altitude levels in GOMETRAN is more 
difficult than the selection of the other computational param- 
eters. The reason for this is that there are effectively two 
degrees of freedom: the number of levels and the choice of 
the altitudes of the levels. In the ultraviolet the selection of 
levels depends on the wavelength, because of rapid changes 
of the optical depth and hence the penetration depth in this 
range. GOMETRAN uses four different altitude grids: three 
in the ultraviolet up to 3 1 5 nm, with a shift of levels from the 
upper to the lower atmosphere toward the longer wavelength 
end; and the fourth grid for the rest of the GOME range with 
a concentration of levels in the lower 35 km of the atmo- 
sphere. In the UV the transitions between different altitude 
grids are at 275 nm, 300 nm, and 315 nm. 

The consequence of a poor selection of the number and 
positions of levels can result in sudden changes in the cal- 
culated radiance, which occur at wavelengths corresponding 
to transitions between the grids, which is demonstrated in 
Figure 4 for the transition occurring at 300 nm. There are 
two different approaches to the determination of the number 
and suitable positions of the atmospheric levels required to 
calculate the radiance accurately. The first is to adopt an al- 
titude grid covering the relevant part of the atmosphere (up 
to about 60 km) with more or less equally spaced levels with 
high vertical resolution and to add lower-resolution levels 
for the remaining height up to about 100 km. A resolution 



Figure 4. Radiance calculated with GOMETRAN before 
optimization of the selected altitude levels. The transition 
between different altitude grids is at 300 nm. Conditions for 
the atmosphere, the viewing geometry, and the computation 
are as in Figure 1 but with four Fourier terms. 


of typically 0.5 km seems to be sufficient and would result 
in about 100 to 150 levels. Because of computational time 
constraints, this is not the preferred long-term solution but 
can be used to establish references for “accurate” values of 
the radiance. 

The second approach to altitude erid selection optimizes, 
for a fixed lower number of levels, the positions of these 
levels in such a way that the accurate reference values cal- 
culated before can be matched by the radiance calculated 
with this lower number of levels. For GOMETRAN it has 
been established that the minimum number of levels should 
be about 50. As the finite difference method requires an odd 
number, 49 has been chosen as the optimum number for the 
altitude levels. The positions of these 49 levels are optimized 
to achieve agreement with an accurate 99-level reference cal- 
culation. 

The choice of altitude levels is dependent on the vertical 
structure of the atmospheric composition. Different scenar- 
ios have been investigated to ensure that a set of optimized 
altitude levels is valid for all climatological conditions. Dif- 
ferences between the radiance calculated with the optimized 
levels and the “reference” radiances (based on 99 levels) are 
displayed in Table 2 for three of the scenarios (Antarctic 
ozone hole condition, northern midlatitude industrial pollu- 
tion in the troposphere, and tropical biomass burning). Also 
presented are the differences between the reference radi- 
ances and those calculated with the original suboptimal lev- 
els (as used in Figure 4). It is evident that the optimization 
has been necessary for a 1 % overall accuracy of the radiance. 
For each of the four wavelength ranges it has been possible 
to find common altitude levels for all scenarios investigated. 
Figure 5 shows for the example of the transition at 300 nm 
the smooth behavior of the radiance achieved after having 
optimized the altitude levels. 

5.5. Radiance at the Line-of-Sight Angle 

The “finite difference method” as described, for example, 
by Lenoble [1985] calculates the radiance only for the zenith 
stream angles used in the Gaussian quadrature of the inte- 
grals in the RTE. Interpolation or extrapolation would then 
be needed to obtain the radiance into the line-of-sight di- 
rection. This process causes considerable errors, especially 
with extrapolation as required with a direct nadir-viewing 
geometry. 

To avoid this error, the line-of-sight angle has been added 
to the regular zenith stream angles with a zero weight The 
result is, because of the zero weight, the maintenance of the 
Gaussian quadrature formula in conjunction with the addi- 
tional calculation of the radiance into difference between 
the reference radiance calculated with 99 altitude levels and 
the radiance calculated with the optimized 49-level grid and 
between the reference radiance and the radiance calculated 
with the original, suboptimal 49-level grid. Conditions for 
the computation are 7 stream zenith angles, 4 Fourier terms 
and 20 Legendre polynomials. 

(Scenario A) Antarctic ozone hole Trace gas profiles 
are based on the MPI database (C. Bruhl , private communi- 
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Table 2. Comparison of Altitude Grids 




Difference (%): 



Difference (%): 




Reference - 



Reference - 




Optimized Grid 



Suboptimal Grid 


X, nm 

Scenario A 

Scenario M 

Scenario T 

Scenario A 

Scenario M 

Scenario T 

274.9 

0.1 

0.3 

0.4 

2.2 

0.2 

0.8 

275.1 

0.2 

0.2 

0.4 

1.5 

2.9 

1.9 

299.7 

0.1 

0.1 

0.8* 

5.0 

5.4 

4.7 

300.1 

0.3 

0.6 

0.1 

0.5 

2.1 

0.7 

314.8 

0.2 

0.4 

0.1 

0.5 

2.1 

0.7 

315.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

789.8 

0.1 

0.1 

0.1 

0.4 

0.4 

0.1 


* TWo additional levels at 3 and 7 km. 

cation (hereinafter referred to as B92)) for 85°S and Octo- 
ber but with a depletion of O3 and NO2 and an enhancement 
of BrO and CIO between 15 and 25 km altitude; maritime 
aerosols between 0 and 2 km, background aerosol types 
above; surface albedo 0.85; solar zenith angle 75°; relative 
azimuth angle between the Sun and the satellite 60°; line-of- 
sight zenith angle 43°. 

(Scenario M) northern midlatitude industrial pollu- 
tion Trace gas profiles are based on the MPI database for 
55°N and January but with elevated levels of O3, NO2, SO2, 
and HCHO in the boundary layer between 0 and 2 km; ur- 
ban aerosols between 0 and 2 km, background aerosol types 
above; surface albedo 0.3; solar zenith angle 75°; relative 
azimuth angle between the Sun and the satellite 0°; line-of- 
sight zenith angle 30°. 

(Scenario T) tropical biomass burning Trace gas pro- 
files are based on the MPI database for 5°N and January 
but with elevated levels of O3, NO2, and HCHO in the tro- 
posphere; rural aerosols between 0 and 2 km, background 
aerosol types above; surface albedo 0.3; solar zenith angle 
30°; relative azimuth angle between the Sun and the satel- 
lite 0°; line-of-sight zenith angle 30°; the direction of the 
line-of-sight angle. The computational time is increased be- 



Figure 5. Radiance calculated with GOMETRAN after op- 
timization of the selected altitude levels, otherwise as in Fig- 
ure 4. 


cause of larger matrix dimensions (by about 30% for seven 
regular stream angles), but the approach is much preferable 
over interpolation/extrapolation. Figure 6 presents a com- 
parison between the radiance at the line-of-sight angle cal- 
culated directly as described above and calculated with an 
interpolation/extrapolation routine based on Legendre poly- 
nomials. It is obvious that errors exceeding the 1% accuracy 
requirement occur if the radiance is not calculated directly. 


To solve the radiative transfer equation, GOMETRAN has 
to solve a system of equations that can be expressed in the 
form of a matrix whose dimensions are roughly equal to 400 
by 400 in case of 50 altitude levels and 8 stream angles. The 
LU decomposition of the matrix performed to obtain this 
solution is the most time-consuming part of the computa- 
tion. GOMETRAN can be operated in two modes: either 
in the monochromatic mode, where the set of equations is 
solved at every wavelength or in the more efficient “quasi- 
polychromatic” mode where a LU decomposition at the sub- 
sequent wavelength is only done when the norm of the new 
matrix has changed by more than a preset factor compared 
to the matrix norm at the previous wavelength. Through ad- 
justment of this factor the frequency of matrix LU decompo- 
sitions can be influenced. The computational speed is high 
for a small number of decompositions but at the expense of 
accuracy and vice versa. 

The O3 Huggins bands between 300 and 350 nm are the 
range with the most rapid changes of the optical proper- 
ties with wavelength. A part of this range was selected 
to optimize the frequency of matrix decomposition. Ref- 
erence values of backscattered radiances were calculated in 
the monochromatic mode with decomposition at each wave- 
length. Differences to the references are shown in Figure 7 
for two conditions, first for decompositions at every 0.2 to 
1 nm, secondly for decompositions every 2 to 10 nm. The 
second case leads to unacceptable errors of up to 3%. The 
first case with the more frequent decompositions represents 
a good compromise, with the accuracy level maintained to 
about 0.5% or lower. 


5.6. Frequency of the LU Decomposition 
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Line-of-Sight Angle [degree] 



Line-of-Sight Angle [degree] 


Figure 6. Direct calculation of the radiance at the line-of- 
sight angle versus interpolation/extrapolation. The normal- 
ized radiance for two wavelengths (top and middle figures) 
is displayed as function of the line-of-sight angle, obtained 
by direct calculation (solid line) and by Legendre interpo- 
lation/extrapolation (dashed line). The seven zenith stream 
angles used as basis in the interpolation/extrapolation pro- 
cess are 13.0°, 29.5°, 45.3°, 60.0°, 72.7°, 82.6°, and 88.5°. 
Also shown is the relative difference between the calculated 
and the interpolated radiances (bottom). Conditions for the 
atmosphere, the viewing geometry, and the computation are 
as in Figure 4 but with the line-of-sight angle being varied. 

6. Comparison with Other Radiative 
Transfer Models 

The plane-parallel GOMETRAN model has been vali- 
dated against other radiative transfer models, especially a- 
gainst DISORT, which is a model based on the “discrete- 
ordinate method.” DISORT in its present plane-parallel form 
is the result of more than 10 years development, the code is 
well documented and validated. The evaluation of DISORT 
is described thoroughly by Stamnes and Swanson [1981], 
Stamnes and Dale [1982], Stamnes and Conklin [1984], 
Stamnes et al. [1988], and Lenoble [1985]. It was shown 


that DISORT is accurate to 1% provided the computational 
parameters are chosen appropriately. 

In DISORT the number of Fourier terms is determined 
automatically by a convergence criterion (which has been 
copied for GOMETRAN). The number of Legendre polyno- 
mials is set internally as twice the number of zenith stream 
angles. Typical numbers of streams (for the half space) in 
the evaluation of DISORT were 8 for moderately anisotropic 
atmospheres, and 16 or 24 when aerosol types are present 
with sharply peaked phase functions. For the aerosol types 
used so far in the GOMETRAN and DISORT test scenarios 
(maritime and tropospheric and stratospheric background 
aerosols), 7 streams were appropriate for both models. Opti- 
mized altitude grids were used for GOMETRAN, while for 
DISORT, optical depth grids with 99 levels were used in the 
ultraviolet and grids with 79 levels in the visible wavelength 
range. 

6.1. Ultraviolet Wavelength Range 

For the comparison between GOMETRAN and DISORT 
in the ultraviolet, radiances were calculated using GOME- 
TRAN in the monochromatic mode. This is necessary to 



305 310 315 320 325 


Wavelength [nm] 

Figure 7. Difference between the GOMETRAN radiance 
calculated with matrix lower-upper (LU) decomposition at 
every GOME wavelength and the radiance obtained with a 
less frequent decomposition: one decomposition every 0.2 
to 1 nm (top), one decomposition every 2 to 10 nm (bottom). 
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avoid unwanted errors caused by infrequent matrix LU de- 
compositions in the Huggins bands between 300 and 350 
nm. The agreement between both forward models is close 
to or better than 0.5% for a set of different viewing geome- 
tries, atmospheric conditions, and wavelength ranges, and it 
is expected that the differences between GOMETRAN and 
DISORT will stay within the 1% range for all possible sce- 
narios. 

Some comparisons were also made with a model devel- 
oped at St. Petersburg University (Af. Birylina and V.V 
Rozanov , private communication, 1992), based on the “dou- 
bling and adding method ” For a northern midlatitude sce- 
nario at a wavelength of 337 nm, a solar zenith angle of 45°, 
a number of line-of-sight and azimuth angles, and surface 
albedos between 0 and 1 the agreement was within 1 %. 

6.2. Visible Wavelength Range 

The radiance calculated by GOMETRAN in the visible 
range shows good agreement with that from DISORT, with 
differences of less than 0.5%. Representative spectral ranges, 
atmospheric and viewing conditions were chosen for the in- 
vestigations. Again, it is anticipated that under all possible 
conditions in the plane-parallel geometry an agreement of 
better than 1 % can be achieved between GOMETRAN and 
DISORT. 

A limited number of comparisons in the visible were made 
between GOMETRAN and two other accurate models, both 
based on the “doubling and adding method.” The first model 
was that from St. Petersburg University, the second one is 
described by Haan et al. [1987] and ( P. Stammes , private 
communication, 1993a). The St. Petersburg comparison 
with GOMETRAN at 650 nm for a northern midlatitude sce- 
nario with a solar zenith angle of 45° and a surface albedo 
of 0.5 produced an agreement of better than 1% for all pos- 
sible line-of-sight and azimuth angles. Comparisons with 
the second doubling and adding model were performed for a 
northern midlatitude scenario and a solar zenith angle of 75°. 
The surface albedo is 0.3, and the line-of-sight and azimuth 
angles are 30° and 60°, respectively. Again, differences of 
less than 1% were observed at 450 and 650 nm. 

7. Conclusion 

A new radiative transfer model, GOMETRAN, has been 
developed for the ultraviolet, visible, and near-infrared wave- 
length range. It has been designed primarily as a tool for the 
data analysis of GOME measurements, but it can also be ap- 
plied to other downward looking spacebome sensors in this 
wavelength range. The model performs a computationally 
efficient calculation of radiances and weighting functions. 
The plane-parallel version of this model, which is applica- 
ble to geographical latitudes between 75°S and 75°N, has 
been presented in this paper. GOMETRAN has been opti- 
mized with respect to the calculation of radiances to achieve 
an accuracy close to 1%. Comparisons with other radiative 
transfer models (based on the discrete-ordinate and on the 
doubling-adding methods) have yielded agreement close to 
or better than 1 % for all parts of the spectral range. 


Future work on GOMETRAN, in addition to the already 
mentioned evaluation of the spherical version and the im- 
plementation of an approximative cloud treatment, will fo- 
cus on two important issues: corrections for the polarization 
insensitivity of the model and the implementation of a pro- 
cedure to take the so-called Ring effect into account. The 
following will briefly discuss these open tasks and possible 
ways for their solution. 

GOMETRAN is an unpolarized model which neglects po- 
larizing properties of the atmosphere in order to limit com- 
putational time requirements. The errors induced by this 
omission have been quantified by a radiative transfer model 
which can be operated both with and without polarization 
taken into account [Haan et ai , 1987], ( P Stammes , private 
communication, 1993b). Typical differences are negligible 
below 290 nm, then rise to a maximum of up to 5% in the 
O 3 Huggins bands, and decrease again toward longer wave- 
lengths with values between 0 and 2 % above 350 nm, de- 
pending on the atmospheric and viewing conditions. To ex- 
ploit the Huggins bands for O 3 profile retrievals and other 
tasks, it will be necessary to correct the radiance calculated 
with GOMETRAN for polarization effects. One possibility 
here is the use of tables with correction factors which have to 
be applied to the unpolarized radiances. These tables would 
have to be created for a variety of atmospheric, surface, and 
viewing scenarios. 

Another point of concern is the Ring effect, which can 
be mainly attributed to inelastic rotational Raman-scattering 
processes in the atmosphere. The most notable feature here 
is the Ailing of the solar Fraunhofer lines [Grainger and 
Ring , 1962]. This effect was shown to be significant in 
the analysis of ground-based spectroscopic measurements, 
where it is usually accounted for by an effective negative 
absorption coefficient [Solomon et al, 1987]. This ignores 
the inelastic nature of the process and is an empirical ap- 
proach based on a limited number of experimental observa- 
tions. The appropriateness of this technique will be stud- 
ied for satellite-based measurements or different approaches 
such as a correction by a single-scattering approximation of 
the Ring spectrum. An accurate treatment is computer in- 
tense and is therefore not considered for GOMETRAN. 
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Abstract — Numerical studies of a new method for the retrieval of ozone profile information 
from nadir-observing satellite measurements in the ultraviolet and visible are presented. The 
method combines information from back scattered radiation in the Hartley band down to the 
O* concentration peak, lower atmospheric information from the temperature structure of the 
Huggins bands, and a constraint on the total column from the Chappuis bands. The Huggins 
bands' temperature structure provides altitude information on the ozone distribution that 
includes clear distinction between stratospheric and tropospheric ozone. Studies presented here 
include dependence of the retrieved 0 3 profiles on surface albedo, tropospheric aerosol, and 
tropospheric 0 3 content for a range of atmospheric conditions. Published by Elsevier Science 
Ltd 


1. INTRODUCTION 

The importance of the stratospheric ozone layer to terrestrial life and the need for precise long-term 
monitoring of the stratospheric ozone climatology are well established. Ozone also plays a 
significant role in the chemistry of the troposphere, where its photolysis leads to the production 
of OH, the most important tropospheric oxidizing agent. 1 Haagen-Smit et al 2 showed that the 
production of tropospheric ozone is dependent on the presence of hydrocarbons, oxides of 
nitrogen, and light. Crutzen 3 elucidated the mechanism for the HO*- and NO*-catalysed 
photochemical generation of ozone in the troposphere. The first experimental evidence that 

tropospheric ozone may be increasing was reported from balloon-sonde observations made 

between 1967 and 1982 at the meteorological station on Hohenpeissenberg in southern Germany. 4 
An analysis of ozone measurements made at the clean air stations of Mauna Loa, Hawaii and Point 
Barrow, Alaska agreed with this finding and indicated an average increase of ozone of 0.5- 1.5% 
per year during the period 1974-1989 with the largest increases occurring during summer. 5 Volz 
and Kley 6 investigated the ozone measurement technique used by Albert-Levy and co-workers at 
Montsouris, near Paris, in the late 19th and early 20th centuries and concluded that the method 
is reliable. 78 These data indicate that the ozone mixing ratio in clean air was about 10 ppbv at the 
end of the last century near Paris. Present day surface measurements are a factor of 2-4 larger. 
As a result of such observations, concern about the impact of tropospheric ozone pollution is 
rapidly growing. To accurately assess the importance of these changes, global measurements of the 
distribution of tropospheric ozone are required. 

Due to the relatively short lifetime and consequent variability of tropospheric ozone, 
measurements made by satellite-borne remote sensing instruments are the most efficient way to 
obtain global information on its distribution. A key difficulty in spacebome monitoring of 

tropospheric ozone is the ability of the measurements and the retrieval approach to distinguish the 
contribution of stratospheric ozone from that of tropospheric ozone in the detected signals. A 
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method to determine the vertical distribution of stratospheric ozone by monitoring back scattered 
radiation in the ultraviolet Hartley and Huggins bands of ozone was first proposed by Singer and 
Wentworth in 1957,’ and has been successfully applied by the BUV instrument on Nimbus 4 and 
the SBUV instrument on Nimbus 7.'° " SBUV obtains stratospheric profile information from the 
differing scattering depths, as well as the total ozone column, by measuring at 12 selected 
wavelength bands in the range 256-340 nm. The SBUV instrument uses relatively wide (1 nm) 
wavelength bands that do not resolve the fine structure of the ozone absorption. Stratospheric 
ozone vertical profiles, with limited geographical coverage, have also been obtained successfully 
using solar occupation by the SAGE I and II instruments. 12 

A method is presented here for remote sensing of ozone profiles, including tropospheric ozone, 
from orbiting spectrometers which measure upwelling atmospheric radiation between 240 and 
800 nm. This study includes simulation of several viewing scenarios, as well as realistic estimation 
of instrument performance derived from studies for the Global Ozone Monitoring Experiment 
(GOME) and the Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY 
(SCIAMACHY), since demonstrated by GOME, to infer the precision of the technique. The 
present work was instrumental in establishing the GOME and SCIAMACHY instrument 
programs, by demonstrating that nadir-looking instruments can measure tropospheric ozone 
globally. GOME is a European Space Agency instrument launched on the European Remote 
Sensing 2 satellite in April 1995. It is a nadir mapping instrument that views the Earth with 
continuous wavelength coverage from 240 to 790 nm at a resolution of 0.2 nm in the ultraviolet 
and 0.4 nm in the visible. 13-15 SCIAMACHY is a German/Dutch/Belgian instrument which is part 
of the payload for the ESA Envisat-1 Polar Platform, currently scheduled for launch in 1999. 
SCIAMACHY is an enhanced version of GOME having, in addition to nadir viewing, limb and 
occultation measuring modes, and extended wavelength coverage to 2.4 /im. 16-18 


2. SPECTROSCOPIC METHOD 

The Huggins bands of ozone have discrete vibrational structure between 300 and 370 nm, with 
features having widths significantly less than 1 nm. This is the sharpest known vibrational structure 
of any electronic band of 0 3 . This structure has a strong temperature dependence due to the onset 
of thermal population of excited vibrational levels in the electronic ground state (see measurements 
in Refs. 19 and 20). It is this temperature-dependent structure that provides a spectral signature 
for tropospheric 0 3 ; its differential character derives from the onset with increasing temperature 
of the weaker vibrational hot-band absorption between the stronger absorption peaks, which are 
due to absorption from 0 3 (000). 

The ozone measurement technique presented here includes the measurement of the Hartley, 
Huggins, and Chappuis bands at moderately high resolution (ca. 0.2 nm in the Hartley and Huggins 
bands, ca. 0.4 nm in the Chappuis bands). In addition to providing more continuous wavelength 
coverage for deriving information on vertical profiles from the scattering depth, it permits the use 
of the temperature structure of the Huggins bands to increase the altitude profile information. In 
particular, since tropospheric ozone is significantly warmer than stratospheric ozone at any altitude 
in the stratosphere where there is significant 0 3 , tropospheric ozone column amounts are cleanly 
separated from the stratospheric column. The Chappuis bands are temperature independent and 
occur at visible wavelengths where light penetrates to the ground under favorable meteorological 
conditions. Thus, their measurement provides improved total column 0 3 information to the 
measurement set. 


3. RADIANCE AND RETRIEVAL CALCULATIONS 
3.1. Forward radiance model 

Radiance calculations are performed for the present studies using the AFGL LOWTRAN7 
radiance code. 21 The version of the code used here includes the correction of all recently discovered 
errata [G. P. Anderson and L. W. Abreu, private communications (1989-1992)]. The only 
modifications made to the code are to increase the number of significant figures for total radiance 
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provided by the code and for ozone concentrations read in by the code. All calculations are 
performed using the provided option for multiple scattering employing the k-distribution method. 22 

Calculations use a 12 layer (13 level) model atmosphere, including 0 3 , H 2 0, and 0 2 , with 
concentrations from the AFGL atmospheric constituent profiles; 23 their model 6 atmosphere, with 
pressure and temperature corresponding to the U.S. Standard Atmosphere, is used. This model 
has a tropopause height of 15 km: subsequent references to tropospheric and stratospheric ozone 
refer to the column from 0 to 15 km and from 15 to 60 km, respectively (the error due to inclusion 
of the small mesospheric column, 50-60 km, with the stratosphere is negligible for present 
purposes). Calculations in this paper are performed for a satellite instrument in a polar orbit with 
a 10:00 a.m. crossing time, looking vertically down, at 45°N, for two viewing scenarios and 
tropospheric 0 3 distributions: (1) vernal equinox, solar zenith angle = 53.5°, ground level 0 3 mixing 
ratio = 40 ppbv (“normal” tropospheric ozone). In this case, the tropospheric ozone column is 63 
Dobson units (DU, 1 DU = 2.687 x 10“ cm -2 ) and the stratospheric column 276 DU; (2) summer 
solstice, solar zenith angle = 48.9°, ground level 0 3 mixing ratio = 120 ppbv (corresponding to a 
boundary layer pollution episode under anticyclonic conditions). In this case, the tropospheric 
ozone column is 76 DU, and the stratospheric column 276 DU. The two ozone distributions are 
shown in Figs l(a-h), along with the retrieved 0 3 results described below. Contributions from 
minor interfering atmospheric species are not included in the present calculations. In the ultraviolet 
these include BrO, H 2 CO, and S0 2 . Measurements of these species simultaneous with ozone 
determinations, implying the appropriate correction of the spectra for analysis of ozone absorption, 
is intended. For example, stratospheric S0 2 , greatly enhanced by the eruption of El Chichon, was 
measured by both the TOMS and SBUV instruments on Nimbus 7. 24 " 26 Clearly, the presence of 
interfering species, even when they are included in the retrieval process, will increase the uncertainty 
in derived ozone over that calculated here. The major interference for the Chappuis bands is from 
N0 2 and H 2 0. These can both be readily corrected for, since the overlaps with the broad Chappuis 
bands are only partial, and since detailed study of the wavelength-dependent absorption by the 
Chappuis bands is not necessary; only their absorption in portions of the bands, cleanly separated 
from interfering species, is necessary to provide column 0 3 measurements. 

The studies include albedos of 0.3 (typical land value) and 0.02 (sea surface minimum) as 
specified in the input to LOWTRAN7. They also include rural and urban tropospheric aerosol 
extinction (IHAZE = 1 and 5 in LOWTRAN7). All studies have been done with the background 
stratospheric aerosol (IVULCN = 0 in LOWTRAN7). Clear sky conditions are assumed. The 
limitations to tropospheric ozone measurements due to cloud coverage are currently being studied 
extensively in the context of the GOME and SCIAMACHY satellite algorithm development. These 
studies include the knowledge of penetration depth of solar radiation into clouds gained from 
measurements in the oxygen A bands. 27 ' 29 

3.2. Instrument modeling 

The instrument model used in the present studies is a simple one developed in the preliminary 
stages of instrument and algorithm design for the GOME and SCIAMACHY projects. The 
instrument is assumed to measure the atmospheric spectrum over the 0 3 regions at 0.2 nm 
resolution from 240 to 350 nm, and at 0.4 nm resolution from 400 to 800 nm (GOME and 
SCIAMACHY both include the 350-400 nm region as well, negligible for 0 3 ; the GOME coverage 
ends at 790 nm, a negligible correction to the coverage in the present study). Spectra calculated 
with LOWTRAN7 are converted to the nanometer scale, appropriate for our measurements with 
dispersive instruments, and are corrected for resolution and for instrument parameters as follows: 
etendue = 5.0 x 10 -5 cm 2 sr; instrument throughput = 0.15 for 2 ^ 400 nm, =0.20 for 2 > 400 nm; 
detector quantum efficiency = 0.5; nadir sampling time = 1.5 sec (the GOME default nadir 
sampling time is 30 sec for 2 < 312 nm and 1.5 sec for 2> 312 nm; for simplicity, and as a 
worst-case scenario, the present study is made using 1.5 sec sampling throughout the GOME 
wavelength range). 

Synthetic noise is added to spectra calculated with the forward model for retrieval study purposes 
using the Box-Muller method for generating random deviates with a Gaussian distribution. 30 Noise 
sources include readout noise (10 3 e' /pixel), dark current noise (10 2 e"/pixel), photon statistics, and 
digitization noise. The digitization noise is calculated assuming 16-bit conversion, with gain factors 
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Fig. I. The results of the fitting study for each of the eight atmospheric scenarios, and for three sets of initial concentrations for each scenario are shown. ( ) 

are the modeled ozone concentrations. Retrieved values include 2 a fitting errors. The atmospheric scenarios are: (a) normal tropospheric ozone, rural haze, albedo 
0.3; (b) normal tropospheric ozone, urban haze, albedo 0.3; (c) normal tropospheric ozone, rural haze, albedo 0.02; (d) normal tropospheric ozone, urban haze, 
albedo 0.02; (e) enhanced tropospheric ozone, rural haze, albedo 0.3; (0 enhanced tropospheric ozone, urban haze, albedo 0.3; (g) enhanced tropospheric ozone, 

rural haze, albedo 0.02; (h) enhanced tropospheric ozone, urban haze, albedo 0.02. 
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set appropriate to GOME measurements under conditions of maximum input. These range from 
1.0 x 10 4 photons per pixel for 268 < A < 312 nm to 9.0 x 10 7 photons per pixel for 
400 ^ X < 800 nm. Figure 2 is an example of the forward calculation, with noise added, 
corresponding to the measurement case with 40 ppbv tropospheric ozone, albedo = 0.3, and the 
rural haze model. 



240 260 280 300 320 340 



600 650 700 750 800 


nm 

Fig. 2. An example of the forward calculation, with noise added, corresponding to the scenario with 
40 ppbv tropospheric ozone, albedo = 0.3, and the rural haze model: (a) the Hartley-Huggins band region. 

Because of the large dynamic range of the spectrum its logarithm is also plotted ( ); (b) the Chappuis 

band region. 
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3.3. Retrieval studies 

Figure 3 shows the ozone absorption cross sections for the Hartley, Huggins, and Chappuis 
bands at a representative temperature for the peak of the stratospheric ozone layer (220 K). These 
are calculated at the GOME/SCIAMACHY spectral resolution. The TOMS and standard (step 




Fig. 3. Ozone absorption cross sections calculated at the GOME/SCIAMACHY spectral resolution, (a) 
The Hartley and Huggins bands at a representative temperature for the peak of the stratospheric ozone 
layer (220 K). The TOMS and standard (step scan mode) SBUV wavelength bands are also shown (TOMS 
also has channels at 360 and 380 nm). The strong variation in the strength for u.v. absorption, together 
with the approximate X 4 dependence of the Rayleigh scattering, provides for discrimination of Oj at 
different altitudes from the measurement of back scattered light at different wavelengths, (b) The Chappuis 

bands. 
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Fig. 4. (a) Temperature dependent cross sections in the O 3 Huggins bands calculated from the 
parameterization in the LOWTRAN7 code over the stratospheric and tropospheric ozone temperature 
range; (b) the same cross sections as Fig. 4(a), but only over the 320-337 nm range, to illustrate details 
of the temperature dependence in the region of primary importance for determination of tropospheric 

ozone. 

scan mode) SBUV wavelength bands are also shown in the figure. The strong variation in the 
strength for u.v. absorption, together with the appropriate X A dependence of the Rayleigh 
scattering, provides for discrimination of 0 3 at different altitudes from the measurement of back 
scattered light at different wavelengths. 

Precise measurements have been made of the temperature dependence of the Huggins bands from 
300 to 340 nm. 19 * 20 The values of the cross sections derived from these measurements currently used 
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in the AFGL LOWTRAN7 radiance code are shown in Fig. 4(a) for the range of temperatures 
relevant to atmospheric ozone measurements. Less precise measurements exist to wavelengths 
longward of 340 nm, in the weaker part of the Huggins bands, and efforts are currently underway 
to provide precise laboratory measurements there as well. Figure 4(b) shows the same cross sections 
as Fig. 4(a), but only over the 320-337 nm range. The portion of the Huggins bands longward of 
320 nm is the region of primary importance for determination of tropospheric ozone. 

Calculated forward spectra, including the instrument model and noise, are generated as described 
in Sec. 3.2. Fitting studies are then performed on these synthetic spectra using the non-linear least 
squares technique. 31 Such an input spectrum is generated, with instrument parameters and noise, 
for each set of input conditions. These include (as detailed in Sec. 3.1) two choices of viewing 
scenario and ozone profile (“normal” and “polluted” troposphere), two choices of albedo (land 
surface average and sea surface minimum), and two choices of tropospheric aerosol distribution 
(rural and urban), for a total of eight input scenarios studied. The non-linear least-squares fitting 
analysis is then performed by allowing the 0 3 concentration at each of the 13 atmospheric levels 
to vary independently in the forward model calculation while iterating to minimize the difference 
between the forward model calculation and the pre-computed forward spectrum plus noise. 
Iteration is continued until the calculated spectrum agrees in the least-squares sense with the input 
spectrum to within a specified x 2 criterion or until all input parameters have stopped varying to 
within preset, small, amounts. Fitting is performed for each of the eight input scenarios, using three 
initial 0 3 distributions. These are: (1) the normal input 0 3 distribution; (2) the normal input 0 3 
distribution multiplied by 0.75, corresponding to 254 DU total ozone; and (3) the normal input 
0 3 distribution multiplied by 1.125, corresponding to 423 DU total ozone. 

Non-linear least-squares fitting of data is not anticipated to be a part of the normal retrieval 
procedure for determination of 0 3 profiles from satellite data; far too much computational time 
is required, even on advanced computers. This type of study has been used during the investigative 
phase of the GOME and SCIAMACHY to determine the maximal information that can be 
obtained from satellite measurements. Algorithm studies now being conducted include, among 
other elements, the use of look-up tables and the optimal estimation method 32 - 33 to design an 
efficient and accurate algorithm for retrieval of atmospheric gas amounts, including 0 3 profiles, 
from the satellite measurements. 


4. RESULTS AND DISCUSSION 

Figures l(a-h) shows the results of the fitting study for each of the eight input scenarios, and 
for the three sets of initial concentrations for each scenario. The concentrations for the input 
spectra and the 2 a fitting errors are included in the figure. Table 1 gives the integrated total ozone 
columns, as well as the tropospheric (i.e., $ 15 km) and stratospheric (> 15 km) columns, and their 
uncertainties, as derived from the fitting analysis. 

The relative fitting uncertainties for tropospheric, stratospheric, and total ozone columns are 
substantially less than those for the individual derived levels. There is substantial correlation 
between ozone derived for adjacent levels, particularly in the troposphere. This is evident from 
examination of the correlation matrices derived in the least-squares fitting procedure (used to 
determine the uncertainties in the column amounts). There is also some suggestion of alternating 
behavior in fitting to adjacent altitude levels in Fig. 1. The fitting to the scenarios with enhanced 
tropospheric ozone shows a tendency to underestimate the ozone at the lowest altitude, 
compensating by placing it higher in the troposphere. The 5 km altitude resolution chosen for 
tropospheric retrievals in this study is probably near to the limit of what is feasible. An exception 
to this may be for ozone in a warm boundary layer where the temperature-induced contrast in the 
Huggins bands will be particularly pronounced. The U.S. Standard Atmosphere temperature and 
pressure profile used in these studies includes a very broad tropopause region (1 1-20 km). When 
atmospheric conditions with steeper temperature gradients in the tropopause region are 
encountered, ozone profiles might be retrievable to higher effective altitude resolution. 

The ability to retrieve ozone profile information depends upon knowledge of the atmospheric 
temperature^ profile. The requirements for this knowledge have not been determined in the present 
study. Temperature knowledge to ± 5 K is estimated to be adequate, but this needs to be quantified 
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Table 1 . Tropospheric stratospheric, and total Oj column retrievals. 



Conditions 

Normal tropospheric ozone t 

Enhanced tropospheric ozonet 

Troposphere 

Stratosphere 

Total 

Troposphere 

Stratosphere 

Total 

“True" profile 

Starting 

Albedo Haze profile} 

62.8 

275.5 

338.3 

76.0 

275.5 

351.5 

0.30 

rural 

L 00 

61.9 + 0.9 

276.1 ±0.6 

337.9 ± 0.3 

75.9 + 0.5 

275.6 ± 0.4 

351.5 + 0.2 

0.30 

rural 

1.25 

59.9+ 1.1 

277.5 ± 0.8 

337.5 + 0.5 

78.1 +0.8 

273.2 + 0.7 

351.2 + 0.4 

0.30 

rural 

0.75 

63.8 ± 0.8 

273.7 ± 0.6 

337.4 ± 0.4 

85.8 ± 0.5 

267.3 ± 0.3 

353.2 ± 0.2 

0.30 

urban 

1.00 

61.7 + 0.6 

276.2 + 0.4 

337.9 + 0.4 

76.0 ± 0.4 

275.5 ± 0.3 

351.5 + 0.3 

0.30 

urban 

1.25 

63.2 + 0.7 

275.2 + 0.5 

338.4 + 0.3 

75.2 + 0.5 

275.8 ± 0.3 

351.0 + 0.3 

0.30 

urban 

0.75 

63.3 ± 0.5 

274.1 ± 0.3 

337.4 ± 0.3 

76.8 ± 0.5 

267.0 + 0.2 

343.8 ± 0.3 

0.02 

rural 

1.00 

62.4 + 0.9 

275.4 ± 0.7 

337.8 + 0.5 

76.5 + 0.5 

275.1 +0.2 

351.6 ±0.3 

0.02 

rural 

1.25 

64.7 ± 1.0 

274.0 + 0.7 

338.7 + 0.5 

83.2 ± 0.6 

267.6 ± 0.3 

350.8 + 0.4 

0.02 

rural 

0.75 

64.3 ± 0.6 

272.8 ± 0.5 

337.1 ±0.4 

82.5 ± 0.6 

268.4 ± 0.4 

350.8 ± 0.5 

0.02 

urban 

1.00 

61.9+ 1.0 

275.5 + 0.7 

337.4 + 0.9 

76.0 + 0.7 

275.5 ± 0.4 

351.5 + 0.6 

0.02 

urban 

1.25 

60.9 ± 2.4 

272.7 ± 0.7 

333.6 + 2.6 

75.5 ± 0.8 

274.6 + 0.5 

350.3 ± 0.7 

0.02 

urban 

0.75 

61.5 ±0.8 

272.6 ± 0.5 

334.1 +0.7 

71.0 ±0.8 

268.2 ± 0.5 

339.3 ± 0.8 


t2 a retrieval errors are given. 

^Multiplier for initial conditions ( x “true” profile with normal tropospheric ozone). 


in further studies. Limited retrievals were also performed using only the Hartley and Huggins band 
region, eliminating the Chappuis bands. In these cases, the uncertainties for tropospheric, 
stratospheric, and total ozone columns were about a factor of 10 higher. An attempt was made 
to retrieve profiles using only the Hartley band but, as expected, there was negligible sensitivity 
to lower atmospheric ozone. 


5. CONCLUSION 

It has been shown that nadir measurements of the ozone Hartley, Huggins, and Chappuis bands 
at sufficient spectral resolution to resolve the structure provides information on the altitude 
distribution of ozone, including tropospheric ozone. Altitude resolution in the troposphere is 
limited to ca. 5 km, except in special conditions. The studies here are limited to clear sky conditions; 
in practice, information could be derived down to the scattering height. These studies contributed 
to the decisions to produce several satellite instruments for atmospheric measurements. The first 
of these, GOME, is now operating successfully in orbit, producing data that should soon confirm 
the level to which this type of retrieval can be carried out practically. 

Table 1 shows that the derived stratospheric and tropospheric ozone columns are cleanly 
separated to quite good precision, even with the caveat given above for strong correlation at the 
5 km altitude resolution scale in the troposphere. Thus, under clear sky conditions, separate 
stratospheric and tropospheric columns are obtainable to high precision to within the limitations 
of the present study. While we do not anticipate that complete global coverage for tropospheric 
ozone will be routinely achieved by the satellite measurements we are undertaking, due chiefly to 
interference in the observations by clouds, we are confident that measurements of its global 
climatology will be greatly enhanced by these satellite measurements. 

This type of measurement and analysis should also substantially improve the determination of 
the climatology of stratospheric ozone, since the interference in measurements by tropospheric 
ozone, including the confusion in analysis caused by the masking of tropospheric ozone by clouds, 
is removed by the clean separation achieved between stratospheric and tropospheric ozone 
determinations. 

The studies for retrieval of ozone and other constituents from the GOME and SCIAMACHY 
instruments are currently in process. Clearly, the work discussed here is preliminary in the sense 
that many further atmospheric and instrumental complications, as well as details of retrieval 
theory, must be considered before a full retrieval procedure can be finalized. Practical application 
may well be limited by systematic effects, such as the uncertainty in radiometric calibration. The 
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results presented here are, however, an encouraging example of the power of utilizing modem array 
detector technology in satellite instruments for atmospheric constituent measurements. These 
efforts should lead to substantial improvement in the determination of the global distributions of 
ozone and other gases. 
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Processing of the RW data to provide a more complete data set is accomplished as follows: 

1. Pressure broadening coefficients for the lines omitted due to interference are replaced by those 
with the same N\ N": R2\Q22 and R25R25 are replaced by R2\R2\ and R25Q26. 

2. 0 2 pressure broadening coefficients measured by RW are multiplied by 1.03, to give air 
broadening. Coefficients for lines measured by RW in the t» = 0«-0 band are also applied to the 
respective lines in the o = l«-l hot band and the u = 0<-0 band of 18 O l6 0. HITRAN92 lines not 
measured by RW are multiplied by 1.03/1.09009 = 0.94488, where 1.09009 is the average ratio of 
the RW measurements to the HITRAN92 coefficients. 

3. The isotopic abundances are taken to be those used in HITRAN92, so that l6 0 2 is 99.52% 
and ,8 0 l6 0 is 0.40%. 

4. Intensity values for the three sub-bands are calculated starting with the 294 K band strength 
of RW, correcting the hot band by the Boltzmann factor of the band origin (1556.385 cm “ l ) and 
the minor isotopic band by the appropriate isotopic ratio. The line intensities are translated to 
296 K by explicitly calculating the rotation-spin partition functions using the term values from 
HITRAN92 and correcting for the RW approximation for the 294 K partition function of the 
major band. 

5. In order to calculate individual line intensities, to supplement the measured ones, it is 
necessary to use the appropriate intensity formulas, which take into account the spin-rotation and 
spin-spin interactions 1 2 8 . For the o = 0<-0 band of 16 0 2 measured by RW the molecular constants 
from Ref. 9 are used, to be consistent, although better constants are now available (Ref. 10 and 
references cited therein). These latter constants are used to calculate the individual line intensities 
for the u = 1«-1 band of 16 0 2 and the u = 0<-0 band of l8 0 16 0. 

The final database derived from these calculation is available from the author via ftp. 

2. TESTING OF LINE SHAPES, INCLUDING LINE NARROWING COEFFICI ENTS 

RW determined that the data as measured under laboratory conditions at high spectral 
resolution ( < 10" 4 cm" 1 ) were fitted better by including narrowing of the line core (Dicke 
narrowing) in the line shape function. Several line shapes were tested, including the Galatry line 
shape M . It was found that all of the tested line shapes which included line narrowing provided 
superior fits to those using the Voigt profile, and that all were equivalent to within experimental 
uncertainty. The Galatry profile was adopted for the bulk of the work because of the relative 
simplicity in implementing it l2 . RW found that all of the lines in the 0 2 A band can be reasonably 
described by a line narrowing coefficient of 0.0145 cm -1 atm -1 (see RW for the exact definition 
of the coefficient in terms of the Galatry profile). An improved method for generating Galatry line 
profiles, employing the fast Fourier transform (FFT), has since been developed by Ouyang and 
Vargese 13 (OV). These authors have distributed a version of their FFT computer code, which was 
used in the present study for investigation of the effect of line narrowing on the atmospheric 0 2 
absorption. In order to use the code, it was necessary to update it and streamline it somewhat. 
The FFT routine used in the code is from the International Mathematical and Statistical Library 
(IMSL); however, calls were to an outdated version of IMSL, which is no longer supported at our 
institution (nor, presumably, at most others). Therefore, the code was updated to IMSL Version 
10. The updated version of the Galatry line profile code is available from the author via ftp. 
Implementation of this code and atmospheric investigation for the 0 2 A band proceeded as follows: 

1. The OV code was compared in the Voigt limit against our standard Voigt subroutine (based 
upon algorithm 363 of the Collected Algorithms from CACM), which has an accuracy of 10 
significant figures, for the case where the Lorentz half-width at half-maximum is 3.0 times the 
Gaussian half-width at 1/e intensity. The maximum disagreement was 0.13%; the maximum 
disagreement relative to the line center intensity was 0.06%. This was accepted as validation that 
the OV code can be properly used to assess the importance of line narrowing on the atmospheric 
0 2 spectrum. 

2. Comparisons were made using the OV code with and without line narrowing for cases 
corresponding to the 0 2 A band at 5 and 10 km altitudes in the atmosphere, using pressure and 
temperature values from the 1976 US Standard Atmosphere l4 . In the 5 km case, the maximum 
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relative error in the line shapes is 1.8%; differences go to < 1% by 2.9 Gaussian 1/e widths from 
line center. In the 10 km case, the maximum relative error is 2.4%; differences go to < 1% by 
2.5 Gaussian 1/e widths from line center. In both cases, the integral of differences over the line 
shapes is zero to within the accuracy of the single precision calculations. In the atmosphere, these 
differences will be in the saturated central parts of the lines for all but the weakest lines. Differences 
greater than 1% are all within 0.02 GOME pixel (for example) so they will not be visible in the 
spectra. Thus, the effect of line narrowing should be totally negligible, and GOME analysis can 
simply use the Voigt profile for satellite-based cloud correction. 




Fig. 1. Template spectra as used in the GOME Initial Cloud Fitting Algorithm, (a) is a sample ground 
template for GOME, i.e., the measured spectrum for light reflected back from ground level 
(1013.25 mbar). (b) is sample template for reflection from a cloud surface at 400 mbar. (c) compares the 
shapes of the two templates to demonstrate the information available in distinguishing the relative level 
of light penetration from satellite spectra. The calculations are made for a solar zenith angle of 60° and 
a satellite viewing angle of 22°, for a total path multiplier of 3.07853. 
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3. SAMPLE 0 2 A BAND CLOUD CORRECTION TEMPLATES FOR GOME 

It has been shown that nadir-looking spectra of the 0 2 A band from space can be used to 
determine cloud-top height and fractional cloud coverage when they are obtained with sufficiently 
high spectral resolution and low noise, and when the database of molecular parameters is 
sufficiently accurate 6 . The GOME satellite instrument 15 initially uses a more approximate 
technique to correct retrieved ozone abundances: Template spectra are prepared for the 0 2 A band 
region for penetration of radiation to the ground and for penetration to standard cloud-top height; 
the effective cloud coverage is retrieved by linear least-squares fitting to determine the best 
coefficients for describing the measurements as a sum of the template spectra. 

For initial comparisons with the GOME satellite data sample ground and 400 mbar cloud 
template spectra are calculated using a 16-layer atmosphere based upon the 1976 U.S. Standard 
Atmosphere 14 and employing Version 5 of the GOME key data for wavelength calibration and 
for the GOME compound hyperbolic slit function (Eq. (1)): 

v= £l + I m 

where y is the normalized spectral response and x is the spectral position in pixel number. The 
current slit function values are a 0 = 0.7334, a { = 0.0756, a 2 = 0.4689, and a 2 - 0.5589. The 
calculations are made for a solar zenith angle of 60° and a viewing angle of 22°, for a total (typical) 
path multiplier of 3.07853, as initially used in the Initial Cloud Fitting Algorithm (ICFA) of 
GOME. Fig. 1(a) is the sample ground template, i.e., the measured spectrum for light reflected 
back from ground level (1013.25 mbar). Fig. 1(b) is the sample template for reflection from a cloud 
surface at 400 mbar. Fig. 1(c) compares the shapes of the two templates to provide some feeling 
for the information available in distinguishing the relative level of light penetration from satellite 
spectra. 
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Ring effect studies: Rayleigh scattering, 
including molecular parameters for rotational 
Raman scattering, and the Fraunhofer spectrum 


Kelly V. Chance and Robert J. D. Spurr 


Improved parameters for the description of Rayleigh scattering in air and for the detailed rotational 
Raman scattering component for scattering by 0 2 and N 2 are presented for the wavelength range 
200-1000 nm. These parameters enable more accurate calculations to be made of bulk molecular 
scattering and of the Ring effect for a variety of atmospheric radiative transfer and constituent retrieval 
applications. A solar reference spectrum with accurate absolute vacuum wavelength calibration, suit- 
able for convolution with the rotational Raman spectrum for Ring effect calculations, has been produced 
at 0.01-nm resolution from several sources. It is convolved with the rotational Raman spectra of 0 2 and 
N 2 to produce an atmospheric Ring effect source spectrum. © 1997 Optical Society of America 


1. Introduction 

The phenomenon that has come to be known as the 
Ring effect was first noted by Grainger and Ring 1 as 
a filling in (broadening and reduction of depth) of 
solar Fraunhofer lines when viewed from the ground 
in scattered sunlight. Various processes have been 
proposed as contributing to the effect, including scat- 
tering with fluorescence from aerosols and from the 
ground . 2 3 However, the predominance of molecular 
scattering as the major cause was established by Kat- 
tawar et aZ ., 4 who analyzed the Ring effect contribu- 
tions from rotational Raman scattering and inelastic 
Rayleighh-Brillouin scattering. The Rayleigh- 
Brillouin component is not of primary importance in 
satellite-based UV-visible backscatter measure- 
ments for which the present study is undertaken, and 
it is not analyzed here. The atmospheric scattering 
situation is nicely defined by Young 6 : “To summa- 
rize: molecular scattering consists of Rayleigh scat- 
tering and vibrational Raman scattering. The 
Rayleigh scattering consists of rotational Raman 
lines and the central Cabannes line. The Cabannes 
line is composed of the Brillouin doublet and the cen- 
tral Gross or Landau-Placzek line. None of the 
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above is completely coherent. The term ‘Rayleigh 
line’ should never be used.” Note that the vibra- 
tional Raman contribution results in lines so widely 
separated from the frequency of the incoming light 
that they are not normally considered part of the 
Ring effect, although in recent applications the Ring 
effect has developed a somewhat broader definition 
that includes substantial interfering structure in ob- 
servations, rather than the initial effect that was lim- 
ited to the broadening of partially resolved lines. 

Understanding the Ring effect has become more 
important in recent years with the increase in ultra- 
violet and visible spectroscopic observations of the 
Earth’s atmosphere from the ground 6 * 7 and from 
satellites . 8-12 To retrieve abundances of trace spe- 
cies from such observations, it is necessary to take 
the Ring effect into account. Methods have been de- 
veloped to do so by pragmatic means: by measuring 
the polarization of scattered sunlight 6 and by model- 
ing the effect directly from molecular scattering pro- 
cesses . 7 - 10 It has been proposed to use the Ring 
effect with selected Fraunhofer lines (in particular, 
the CanH and K lines) to determine cloud parame- 
ters in conjunction with satellite-based observations 
of O 3. 1 1 Joiner et aL 10 also determined that the con- 
tribution to the Ring effect from the Rayleigh- 
Brillouin scattering process is negligible for most 
geometries used in satellite observations. 

The use of visible bands of 0 2 , in particular the 
762-nm A band, for determination of cloud parame- 
ters is being developed for proposed and existing sat- 
ellite instruments that monitor the atmosphere with 
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Table 1. Relative Rayleigh and Raman Scattering Intensities* 


V Polarization in 

H Polarization in 

Sum (Natural Light in) 

Rayleigh-Brillouin 



V C V = 180 + 4s 

H Cy = 3S 

°Cy = 180 + 7s 

p' 

II 

CO 

m 

h C h = 3e + (180 + e)cos 2 0 

°C H = 6s + (180 + s)cos 2 0 

V C 0 = 180 + 7s 
Raman 

H C 0 = 6e 4- (180 + e)cos 2 0 

°C 0 = (180 + 13s) + (180 + s)cos 2 8 
p 0 c = 6s/(180 + 7s) 

V W V = 12s 

X = 9s 

°W V = 21s 

II 

SP 

h W h = 9s + 3s cos 2 0 

°W„ = 18s + 3s cos 2 0 

V W 0 = 21s 
Sum 

H W 0 = 18s + 3s cos 2 9 

0 Wq = 39s + 3s cos 2 8 
Po W = 6/7 

v Ty = 180 + 16s 

H T V = 12s 

°T V = 180 + 28s 

V T„ = 12s 

h T h = 12s + (180 + 4e)cos 2 9 

°T„ = 24s + (180 + 4 s)cos 2 6 

V T 0 = 180 + 28s 

H T 0 = 24s + (180 + 4s)cos 2 e 

°T 0 = (180 + 52s) + (180 + 4 s)cos 2 6 
Po T = 6s/ (45 + 7s) 


“Mostly from Kattawar et al . , 1981. 


emphasis on tropospheric measurements . 13 This in- 
cludes the European Space Agency's global ozone 
monitoring experiment (GOME) and the upcoming 
scanning imaging absorption spectrometer for atmo- 
spheric chartography (SCIAMACHY). The Ring ef- 
fect has a substantial influence on such observations; 
work is under way in our institution to refine and 
model the effects. In this case, the effect is due to 
inelastic scattering in molecular absorption lines 
themselves rather than in the solar Fraunhofer lines. 
S imil ar effects have been noted for the detailed re- 
trieval of trace species including 0 3 and N0 2 . 6,7 12 

This study is part of an ongoing effort to quantify 
the Ring effect for atmospheric radiative transfer 
modeling, with application to satellite- and ground- 
based measurements, and to apply it to particular 
cases such as detailed absorption in the visible 0 2 
bands, the retrieval of tropospheric 0 3 from UV mea- 
surements, and retrieval of trace photochemically ac- 
tive gases. Much of the previous modeling 
work 7 - 10 * 14 has relied on the development of molecular 
parameters for N 2 and 0 2 by Penney et a/. 15 ; only one 
study 10 has used the updated dynamic polarizability 
anisotropies as developed by Bates . 16 Previous re- 
search has largely ignored the complication of the 
rotational Raman spectrum of 0 2 caused by the elec- 
tronic spin angular momentum in the ground 
state and the issue of pressure broadening of the 
rotational Raman lines. In this publication, we up- 
date the molecular parameters and the scattering 
with respect to the solar Fraunhofer spectrum using 
the best currently available laboratory and field data 
and theoretical studies of which we are aware. This 
provides an updated expression for Rayleigh scatter- 
ing by air, expressions for the wavelength-dependent 
polarizability anisotropies of 0 2 and N 2 , accurate 
Placzek-Teller coefficients (the state-dependent fac- 
tors in the line intensities) for 0 2 rotational Rama n 
lines (the Placzek-Teller coefficients of Penney et al. 15 
for N 2 are correct as given), a tentative set of 
pressure-broadening coefficients for the 0 2 and N 2 
rotational Raman lines, a solar reference spectrum 


for convolution with calculated Ring cross sections, 
and a convolved Fraunhofer-rotational Raman source 
spectrum for fitting of atmospheric spectra. The ta- 
bles and spectra are not included here because of size 
limitations, but they are available from the authors. 

2. Rayleigh Scattering 

To e xamin e the detailed Rayleigh and rotational Ra- 
man scattering properties, including their relative 
intensities and the scattering phase functions, we 
begin with Table I of Kattawar et al., 4 which de- 
scribes the relative intensities and angular behaviors 
for the Rayleigh-Brillouin and rotational compo- 
nents and their sum. They are given for various 
input polarizations, including unpolarized light (the 
predominant contribution for many atmospheric ob- 
servations for which single scattering is the major 
contributor to the Ring effect). For unpolarized 
light, the depolarization ratios (defined in each case 
as the ratio of the horizontally polarized component 
to the vertically polarized component at a 90° scat- 
tering angle) can be determined directly. Table I of 
Kattawar et al 4 is reproduced here as Table 1 , with 
the addition of depolarization ratios for three cases: 
Rayleigh-Brillouin (the central Cabannes component 
C); rotational Raman (the wings W); and the sum of 
the two (T for total). The phase functions for scat- 
tering can also be derived for each case. They are 
(normalized over solid angle to 1 ) 

(180 + 13e) + (180 + e)cos 2 8l 


* o C = 


160tt 


V = ^ (13 + cos 2 e), 

(45 + 13e) + (45 + €)cos 2 9l 
9 + 2e J ’ 

where c = ( 7 /a) 2 , 7 is the anisotropy of the polariz- 
ability, and a is the average polarizability. In each 
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case, the phase function is given in terms of the re- 
spective depolarization ratio p 0 * (X = C,W, T) by 


3 [ (1 + p 0 x ) + (1 ~ Pq*)cos 2 9 
8-rr 2 + p 0 x 


(2) 


The Rayleigh scattering cross section at standard 
temperature (273.15 K) and pressure (1 atm) is given 
by 


Qr “ 


32t r 3 (n - 1 fF K 
3N Q 2 \ 4 


(3) 


where n is the index of refraction, F K is the King 
correction factor, 17 N 0 is Loschmidt’s number 
(2.686763 x 10 19 cm -3 ), and \ is the wavelength. 
The King correction factor is given by 


King factor and depolarization ratio for any choice of 
wavelength within the 0.2-1. 0-|xm range can be accu- 
rately determined from the previous two equations. 

B. Polarizability Anisotropies 

The previous expansions are for air, including stan- 
dard amounts of Ar and C0 2 . For rotational Raman 
cross sections, we limit the calculations to 0 2 and N 2 . 
Bates gave segmented representations for the indices 
of refraction of 0 2 and N 2 versus wavelength and 
expansions for the respective King correction factors. 
Tables giving the King correction factors, values of 
Ve, index of refraction a, and y for 0 2 and N 2 have 
been derived in this study. Values of y 0z from these 
data are fitted to 0.6% rms over this wavelength 
range (1% worst case) by 


1 + 2|-^| = 

oct 


6 + 3p 0 7 


2e 


~ 7 T = 1 + IT ' 

6 - 7p 0 9 


The average polarizability can be determined from 
the relation 


91 45.9364 

7o 2 x 10 24 = 0.07149 + 48 2?16 _ ~i • (8) 

The 7 jv 2 values are fitted to 0.02% rms (0.04% worst 
case) by 


Id| 2 = 


(ft ~ l) 2 
WN 0 2 ' 


(5) 


1 n 2 X 10 25 = -6.01466 + 


2385.57 
186.099 - cr 2 ' 


(9) 


3. Molecular Parameters 

A. Rayleigh Scattering Cross Sections 

The major source of improved data for the indices of 
refraction of 0 2 , N 2 , and air versus wavelength, the 
King correction factors, and the anisotropies of the 
polarizabilities is Bates. 16 He presented a compre- 
hensive review of both measurements and theoretical 
calculations and derived a data set that was demon- 
strated to be better than 1% for all the data and 
parameterizations utilized in the this study. 

Bates's Table 1 gives refractive indices, Rayleigh- 
scattering cross sections, and King correction factors 
versus wavelength for air from 200 to 1000 nm. We 
found that the following Edlen-type expression fits the 
index of refraction data to better than 0.1% for all 
values: 


(n air - 1) x 10 4 = 0.7041 + 


315.90 
157.39 - cr 2 


8.4127 

+ 50.429 - a 2 ’ 


(6) 


where cr(|xm -1 ) = 1/X (*xm). The Rayleigh cross sec- 
tions are reproduced to 0.3% rms (worst case, 0.5%) 
over the full 0.2-l-|xm range by the expression 


Q R x 10 24 (cm 2 ) 

3.9993 x 10'V 

“ 1 - 1.069 x 10" V - 6.681 x lO'V * 


(7) 


Nicolet 18 gave an expression for fitting the data of 
Bates to 0.5% over the 0.2-0.55-pm range; the current 
expression is slightly more accurate and extends over 
a larger wavelength range (also see Bucholtz 19 ). The 


These equations should be considered merely as phe- 
nomenological fits over a particular wavelength 
range, rather than trying to attach physical impor- 
tance to, for example, the negative value of the con- 
stant term in the N 2 equation. They should not be 
extrapolated to outside the 0.2-1. 0-^m range. 

C. Basic Spectroscopy 

The ground state of 0 2 is 3 2~; it has significant elec- 
tronic structure in both its magnetic dipole rotational 
spectrum and its rotational Raman spectrum. Ro- 
tational Raman spectra of both 0 2 and N 2 have pre- 
viously been approximated by simple expansions in 
the lowest rotational parameters (for positions) and 
by using T/c<^B Q for the rotational partition functions. 
Precise data are now available for the term energies, 
allowing line positions and Boltzmann factors to be 
calculated accurately and rapidly. The use of term 
values from the current hitran listing 20 allows for 
calculations of line positions to 0.0001-cm -1 accuracy 
and for highly accurate calculation of statistical par- 
titioning. 

D. Placzek-Teller Coefficients 

Cross sections for rotational Raman scattering are 
given by 

Q n , n w ( cm 2 ) = y%cMN, J, N', J'), (10) 

where f N is the fractional population in the initial 
state. The quantum state-dependent factors in the 
cross sections CpjiN, J, N\J r ) are commonly known as 
Placzek-Teller coefficients from the initial deriva- 
tions. 21 These factors are given for a molecule in a 2 
electronic state with one electronic spin angular mo- 
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mentum in the Hund’s case b coupling scheme (elec- 
tron spin coupled to rotational angular momentum) by 

CpxW J, AT, J ’) = (2 N+ 1)(2 N' + 1)(2 J’ + 1) 


l N 

L 

N'\ 2 jN 

L 

API 

\0 

0 

o) W' 

S 

A 


(ID 

where J is the total angular momentum, S is the 
electronic spin angular momentum (1 for 0 2 and 0 for 
N 2 ), and L is the component of the 2nd rank polariz- 
ability tensor responsible for the interaction. L - 0 
for the isotropic part of the polarizability (the Ca- 
bannes component), and L = 2 for the anisotropic 
part of the polarizability (Raman scattering). The 
standard definitions for 3 -j and 6^‘ coefficients are 
used. Two near equivalents to this equation give 
the formula for line strengths. These include the 
initial state degeneracy. 22 * 23 For S = 0 (i.e., N 2 ), Eq. 
(11) reduces to the result of Eq. (7) in Ref. 15: 


— *«/ "H 2) 
Cpr(J -»J) 


3 (J + l)(e/+2) 
2(2J + 1)(2 J + 3) * 

J(J+ 1) 

(2 J - 1)(2 J + 3) ’ 


Cpj^ J — 2) — 


3J(J- 1) 

2(2J + 1)(2J - 1) * 


(12) 


This derivation gives correct Placzek-Teller coeffi- 
cients for N 2 but approximates those for 0 2 by treating 
it as a pure Hund’s case b molecule. At low N and J , 
the departures from the pure coupling case that are 
due to the electron spin-rotation interaction are 
enough to affect the spectrum significantly. 22 * 24 In 
the classic study of the 0 2 ground state, Tinkham and 
Strandberg 25 include the correct eigenvectors of the 
secular determinant for levels to as high as J = 26 in 
the case b basis set (their Table V and Eq. 54). Above 
this level, the molecule is described by case b behavior 
to a high degree of accuracy. For J = 0 and all odd J 
levels, the case b description is exact. For even J 
levels, J =£ 0, the transformation from case b basis 
functions to eigenfunctions is given by 


= dj§j-u + bj$j+ ij. (13) 

The b j and dj values from Ref. 25 are used here to 
calculate correct Placzek-Teller coefficients for val- 
ues to as high as J = 10, above which the corrections 
become completely negligible. An almost identical 
result was determined by Altmann et a/., 26 who cal- 
culated eigenvectors for the secular determinant and 
molecular parameters given in the slightly earlier 
work of Mizushima and Hill. 27 The database calcu- 
lated with the corrected eigenvectors is prepared with 
an intensity cutoff to include all rotational Raman 
lines with intensities at 296 K within 0.1% of the 
strongest line. Because of the mixing of states, this 




Fig. 1. Rotational Raman cross sections for (a) N 2 and (b) 0 2 
calculated for a temperature of 250 K and an excitation wavelength 
of 440 nm, appropriate to investigations of the effect of the Ring 
effect on retrievals of atmospheric NO z concentrations. The cross 
sections in this figure do not include the relative abundances of N 2 
and 0 2 in air. 


now includes two AAT = 4 transitions. Figure 1 gives 
rotational Raman cross sections for N 2 and 0 2 calcu- 
lated for a temperature of 250 K and an excitation 
wavelength of 440 nm, appropriate to investigations 
of the effect of the Ring effect on retrievals of atmo- 
spheric N0 2 concentrations. 

E. Pressure-Broadening Coefficients 

The capital letter T is used here to denote the half- 
width at half-maximum pressure-broadening coeffi- 
cient to distinguish it from the y used for anisotropy 
of polarizability. The best existing measurements of 
pressure broadening for the rotational Raman lines 
are from Jammu et al. 28 who noted some evidence 
that the unresolved Q branches of vibrational Raman 
bands seem to broaden less than lines caused by or- 
dinary dipole transitions but that the rotational Ra- 
man lines broaden comparably to dipole transitions. 
They present self-broadening measurements for N 2 , 
C0 2 , and CO, as well as He and Ar broadening of 
N 2 , 0 2 . All the measurements were made at room 
temperature. Given the experimental conditions 
(high pressures, modest spectral resolution) and the 
lack of air-broadening measurements, temperature 
dependences, and a tabulation of the J-dependent N 2 
broadening coefficients, it was decided not to use 
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Fig. 2. (a) Sample of the Fraunhofer reference spectrum determined in this study, (b) The same Fraunhofer spectrum, convolved with 

a Gaussian slit function having a half-width at 1/e intensity of 0.212 nm, appropriate to the GOME instrument in this wavelength range, 
(c) The Fraunhofer spectrum of (b), convolved with the rotational Raman cross sections to create a Ring effect source spectrum. The units 
are those of cross section x photons, corresponding to the Ring effect scattering source per air molecule at 250 K (d) The ratio of the Ring 
effect source spectrum given in (c) to the Fraunhofer spectrum of (b), with a cubic polynomial subtracted off. This spectrum closely 
corresponds to the Ring effect spectrum used in previous studies for atmospheric spectrum correction, as discussed in the text. 


these measurements as the basis for broadening in 
the present data set. 

For 0 2 , the hitran 92 values corresponding to ro- 
tational transitions of the same AN are initially 
adopted. 20 When multiple corresponding transi- 
tions exist, the pressure-broadening coefficients are 
averaged. For the two A N = 4 transitions, the av- 
erages of the values for lines connecting the upper 
and lower states are taken. The resulting pressure- 
broadening coefficients are multiplied by 1.185, which 
is the average result for the ratio of measured air 
pressure-broadening coefficients for 0 2 magnetic di- 
pole rotational transitions to those given in hitran 
(see Ref. 29 for an explanation of this correction). 
For N 2 , pressure-broadening values for the corre- 
sponding quadrupole lines of the vibrational funda- 
mental are adopted. Values determined for the 
pressure-broadening coefficients of both 0 2 and N 2 
are at 296 K. The temperature dependence should 
be calculated using the recommended hitran 92 coef- 
ficient of n = 0.5 for N 2 lines 20 and a value of n — 0.72 
for 0 2 lines, 29 through the relation 



4. Solar Reference and Ring Source Spectra 

A solar reference spectrum for the range 230-800 
nm, at 0.01 nm resolution, has been determined by 
combining ground-based measurements 30 and bal- 
loon measurements. 31 Both spectra were converted 
to vacuum wavelengths. The balloon data were 
recalibrated in wavelength using 20 selected atomic 
reference lines. The spectra were resampled at even 
0.01-nm increments, employing a triangular filter of 
0.01-nm half-width at half-maximum, and linearly 
merged over the 300-305-nm wavelength range. 
The result is a reference spectrum in vacuum wave- 
lengths accurate to better than 0.001 nm above 305 
nm and 0.002 nm below 300 nm. This spectrum was 
created specifically for GOME applications, in partic- 
ular wavelength calibration studies and calculation 
of Ring effect contributions to GOME measurements. 
A sample of this Fraunhofer reference spectrum is 
given in Fig. 2(a); in Fig. 2(b), it has been convolved 
with a Gaussian slit function having a half-width at 
1/e intensity of 0.212 nm, appropriate to the GOME 
instrument in this wavelength range. The solar 
spectrum, with GOME slit function, has been con- 
volved with the rotational Raman cross sections de- 
scribed here to create a Ring effect source spectrum 
for fitting of GOME data. The same wavelength 
portion of this Ring effect source spectrum is shown in 
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Fig. 2(c). The units are those of cross sections times 
photons, corresponding to the Ring effect scattering 
source per air molecule at 250 K. In Fig. 2(d), the 
ratio of the Ring effect source spectrum to the Fraun- 
hofer spectrum has been taken and a cubic polyno- 
mial, fitted to the region presented in the figure, has 
been subtracted off, for comparison with differential 
Ring effect spectra determined in previous studies. 
Atmospheric Ring spectra have been determined by 
several groups from polarized spectroscopic measure- 
ments of the zenith sky at high solar zenith angles. 
Measurements at polarizations parallel and perpen- 
dicular to the single-scattering plane are combined in 
similar formulations to obtain effective Ring effect 
cross sections. 6 ’ 742 For a single-scattering Rayleigh 
atmosphere, this procedure could yield the rotational 
Raman cross sections directly. In practice, because 
of Mie and multiple Rayleigh scattering, what is de- 
termined is the ratio of the rotational Raman spec- 
trum to the solar Fraunhofer spectrum. This can be 
confirmed by comparing Fig. 2(d) with the results of 
the above studies. Atmospheric radiative transfer 
calculations of Ring spectra 712 also determine this 
quantity. It is distinct from the filling in calculated 
by Joiner et aZ., 10 who include both the rotational 
Raman source term and the Cabannes loss term in 
their calculations. The use of the ratio of rotational 
Raman to Fraunhofer has proved useful in fitting 
atmospheric spectra for minor constituent abundanc- 
es; it is a good approximation to the quantity needed 
for Ring effect correction in some fitting techniques 
and measurement geometries. However, its use has 
limi tations because of the distortion of the rotational 
Raman spectrum and its constant ratio to the Fraun- 
hofer spectrum. The rotational Raman spectrum 
and Fraunhofer spectrum determined in this study 
should be more generally applicable to fitting of at- 
mospheric spectra obtained from various measure- 
ment geometries, including satellite-based 
measurements. 

5. Conclusions 

The complete Rayleigh with rotational Raman- 
scattering database as described in the previous sec- 
tions is available from the authors. It includes N 2 
and 0 2 term values and statistical partitioning (ta- 
bles including the quantum numbers and term ener- 
gies for 0 2 and N 2 up to states allowing for 
partitioning to better than 0.01% accuracy and Boltz- 
mann factors with nuclear spin degeneracies^); ta- 
bles giving the King correction factors, values of Ve, 
where € = (y/a) 2 , index of refraction a, and y for 0 2 
and N 2 ; Placzek-Teller coefficients, with examples of 
cross sections; pressure-broadening coefficients; the 
solar reference spectrum derived above; and the con- 
volved solar and rotational Raman spectrum (Ring 
effect source function cross section calculations). 
This database is a summary of the best available 
relevant data that were found in the literature and 
the calculations performed in this study. One item 
not included here is the additional broadening that is 
due to Rayleigh-BriUouin scattering. 4 For the rota- 


tional Raman lines, this will provide an extra source 
of broadening, although the extent of the broadening 
is within the uncertainties in the pressure broaden- 
ing for scattering in the troposphere. The effect of 
Rayleigh-BriUouin scattering for the fiUing in of the 
central Cabannes line for narrow Fraunhofer lines 
might be significant for some satellite measurement 
conditions and wifi be a topic for future investiga- 
tions. 
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